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ABSTRACT 


(r) 

A  study  was  made  of  radiation-induced  dimensional  changes  in  Pyrex  and  Hoya 
SD-2^  glasses.  These  glasses  are  used  as  substrates  for  MEMS  devices  employing  silicon 
oscillating  beams,  and  changes  in  substrate  dimension  can  compromise  device  performance. 
Silicon  MEMS  strain  gauges  mounted  on  glass  substrates  were  exposed  to  gamma  doses  up 
to  385  Mrad.  The  device-substrate  differential  strain  thus  measured  corresponds  to  an 
expansion  in  the  glass  following  a  linear  trend  with  dose  of  about  5  (±4)  X  10  ^  V tad  and  4.5 
(±0.2)  X  /rad  for  Pyrex  and  Hoya  SD-2,  respectively.  Separate  glass  samples  were 
irradiated  with  a  neutron  fluence  composed  primarily  of  thermal  neutrons,  and  a  floatation 
technique  was  employed  to  measure  the  resulting  density  changes  in  the  glasses  alone. 
Transport  of  Ions  in  Matter  (TRIM)  calculations  were  performed  to  relate  fast  (~lMeV) 
neutron  atomic  displacement  damage  to  diat  from  boron  thermal  neutron  capture  events, 
and  measured  density  changes  in  the  glass  samples  were  thus  proportionally  attributed  to 
thermal  and  fast  neutron  fluences.  Maximum  neutron  fluences  received  by  glass  samples 
were  estimated  at  4  X  10*^’  n/cxn  thermal  and  6  X  n/cra  fast.  The  glasses  irradiated 
with  thermal  neutrons  were  found  to  exliibit  Unear  compaction  at  a  rate  of -2.8  X  10  per 
n/crn  for  Pyrex  and  -1.0  X  10  per  n/cm  for  Hoya  SD-2.  For  fast  neutron  fluence,  strain 
rate  (per  n/cxn)  was  also  Unear:  -6.1  X  10^^  for  Pyrex  and  -7.9  X  10^^  for  Hoya  SD-2.  The 
gamma  irradiation  strain  gauge  results  cannot  be  explained  by  the  radiation-induced 
compacdon  previously  reported  for  Pyrex,  and  are  thought  to  result  from  either  changes  in 
glass  radiadon  response  caused  by  the  anodic  bonding  process  used  to  make  the  MEMS 
strain  gauges  or  from  radiadon-assisted  creep  of  the  strain  gauge  anchors.  The  neutron 
irradiadon  results  for  Pyrex  show  good  agreement  with  those  of  previous  studies. 
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ABSTRACT 

A  study  was  made  of  radiation-induced  dimensional  changes  in  Pyrex  and  Hoya 
SD-2®  glasses.  These  glasses  are  used  as  substrates  for  MEMS  devices  employing  silicon 
oscillating  beams,  and  changes  in  substrate  dimension  can  compromise  device  performance. 
Silicon  MEMS  strain  gauges  mounted  on  glass  substrates  were  exposed  to  gamma  doses  up 
to  385  Mrad.  The  device-substrate  differential  strain  thus  measured  corresponds  to  an 
expansion  in  the  glass  following  a  linear  trend  with  dose  of  about  5  (±4)  X  10  ’^  / rad  and  4.5 
(±0.2)  X  10  /rad  for  Pyrex  and  Hoya  SD-2,  respectively.  Separate  glass  samples  were 
irradiated  with  a  neutron  fluence  composed  primarily  of  thermal  neutrons,  and  a  floatation 
technique  was  employed  to  measure  the  resulting  density  changes  in  the  glasses  alone. 
Transport  of  Ions  in  Matter  (TRIM)  calculations  were  performed  to  relate  fast  (~lMeV) 
neutron  atomic  displacement  damage  to  that  from  boron  thermal  neutron  capture  events, 
and  measured  density  changes  in  the  glass  samples  were  thus  proportionally  attributed  to 
thermal  and  fast  neutron  fluences.  Maximum  neutron  fluences  received  by  glass  samples 
were  estimated  at  4  X  n/ctn  thermal  and  6  X  lO'"*  n/cm^  fast.  The  glasses  irradiated 
with  thermal  neutrons  were  found  to  exhibit  linear  compaction  at  a  rate  of  -2.8  X  10'^'’  per 
n/cm^  for  Pyrex  and  -1.0  X  10'^'  per  n/crn  for  Hoya  SD-2.  For  fast  neutron  fluence,  strain 
rate  (per  n/ cm^)  was  also  linear:  -6.1  X  10'^'  for  Pyrex  and  -7.9  X  10  for  Hoya  SD-2. 
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ABSTRACT 

A  study  was  made  of  radiation-induced  dimensional  changes  in  Pyrex  and  Hoya 
SD-2®  glasses.  These  glasses  are  used  as  substrates  for  MEMS  devices  employing  silicon 
oscillating  beams,  and  changes  in  substrate  dimension  can  compromise  device  performance. 
Silicon  MEMS  strain  gauges  mounted  on  glass  substrates  were  exposed  to  gamma  doses  up 
to  385  Mrad.  The  device-substrate  differential  strain  thus  measured  corresponds  to  an 
expansion  in  the  glass  following  a  linear  trend  with  dose  of  about  5  (±4)  X  10  / rad  and  4.5 

(±0.2)  X  10  ’'*  /rad  for  Pyrex  and  Hoya  SD-2,  respectively.  Separate  glass  samples  were 
irradiated  with  a  neutron  fluence  composed  primarily  of  thermal  neutrons,  and  a  floatation 
technique  was  employed  to  measure  the  resulting  density  changes  in  the  glasses  alone. 
Transport  of  Ions  in  Matter  (TRIM)  calculations  were  performed  to  relate  fast  (~lMeV) 
neutron  atomic  displacement  damage  to  that  from  boron  thermal  neutron  capture  events, 
and  measured  density  changes  in  the  glass  samples  were  thus  proportionally  attributed  to 
thermal  and  fast  neutron  fluences.  Maximum  neutron  fluences  received  by  glass  samples 
were  estimated  at  4  X  lO'*^  n/cm^  thermal  and  6  X  10’"^  n/cm^  fast.  The  glasses  irradiated 
with  thermal  neutrons  were  found  to  exhibit  linear  compaction  at  a  rate  of  -  2.8  X  10  ""  per 
n/cm"  for  Pyrex  and  -1.0  X  lO  "’  per  n/cm"  for  Hoya  SD-2.  For  fast  neutron  fluence,  strain 
rate  (per  n/cm")  was  also  linear:  -6.1  X  10  "'  for  Pyrex  and  -7.9  X  10""  for  Hoya  SD-2. 
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ABSTRACT 

A  study  was  made  of  radiation-induced  dimensional  changes  in  Pyrex  and  Hoya 
SD-2®  glasses.  These  glasses  are  used  as  substrates  for  MEMS  devices  employing  silicon 
oscillating  beams,  and  changes  in  substrate  dimension  can  compromise  device  performance. 
Silicon  MEMS  strain  gauges  mounted  on  glass  substrates  were  exposed  to  gamma  doses  up 
to  385  Mrad.  The  device-substrate  differential  strain  thus  measured  corresponds  to  an 
expansion  in  the  glass  following  a  linear  trend  with  dose  of  about  5  (±4)  X  10'^^  /rad  and  4.5 
(±0.2)  X  10  /rad  for  Pyrex  and  Hoya  SD-2,  respectively.  Separate  glass  samples  were 
irradiated  with  a  neutron  fluence  composed  primarily  of  thermal  neutrons,  and  a  floatation 
technique  was  employed  to  measure  the  resulting  density  changes  in  the  glasses  alone. 
Transport  of  Ions  in  Matter  (TRIM)  calculations  were  performed  to  relate  fast  (^^IMeV) 
neutron  atomic  displacement  damage  to  that  from  boron  thermal  neutron  capture  events, 
and  measured  density  changes  in  the  glass  samples  were  thus  proportionally  attributed  to 
thermal  and  fast  neutron  fluences.  Maximum  neutron  fluences  received  by  glass  samples 
were  estimated  at  4  X  10^^’  n/cvn  thermal  and  6  X  10^"^  n/crn  fast.  The  glasses  irradiated 
with  thermal  neutrons  were  found  to  exhibit  linear  compaction  at  a  rate  of  -2.8  X  10'^^  per 
n/crn  for  Pyrex  and  -1.0  X  10  ^^  per  n/crn  for  Hoya  SD-2.  For  fast  neutron  fluence,  strain 
rate  (per  n/ cm^)  was  also  linear:  -6.1  X  10’^^  for  Pyrex  and  -7.9  X  10  ^  for  Hoya  SD-2.  The 
gamma  irradiation  strain  gauge  results  cannot  be  explained  by  the  radiation-induced 
compaction  previously  reported  for  Pyrex,  and  are  thought  to  result  from  either  changes  in 
glass  radiation  response  caused  by  the  anodic  bonding  process  used  to  make  the  MEMS 
strain  gauges  or  from  radiation-assisted  creep  of  tlie  strain  gauge  anchors.  The  neutron 
irradiation  results  for  Pyrex  show  good  agreement  with  those  of  previous  studies. 

Technical  Supervisor:  Dr.  Jeffrey  Borenstein 

Title:  Principal  Member  of  the  Technical  Staff 


Thesis  Advisor: 
Title: 


Professor  Linn  W.  Hobbs 
Professor  of  Materials  Science 


ACKNOWLEDGEMENTS 

August  4,  2000 


Thanks  go  to  the  people  in  the  United  States  Air  Force  who  decided  to  give  me  an 
opportunity  for  schooling.  Upon  my  return  to  a  "real"  Air  Force  job,  I  hope  my  job 
performance  will  prove  them  right  in  supporting  higher  education  for  Air  Force  people. 
Thanks  also  go  to  the  education  office  folks  (|ohn  Sweeney,  George  Schmidt,  Loretta 
Mitrano,  and  Michael  Villani)  for  running  the  Draper  Fellow  program,  without  which  I 
wouldn't  be  here  at  all. 

This  thesis  is,  of  course,  the  result  of  efforts  by  many  people  not  included  on  the  tide 
page.  |eff  Borenstein  and  Linn  Hobbs  deserr^e  credit  for  guiding  me  during  the  whole 
process.  It  was  they  who  came  up  with  the  idea  of  using  the  pointer  beam  strain  gauges  in 
the  radiation  experiments.  Prof  Hobbs  furnished  me  with  his  analysis  of  the  existing 
literature,  and  suggested  the  TRIM  calculations.  Dr.  Borenstein  arranged  for  the  fabrication 
of  the  pointer  beams  and  signed  purchase  requests  in  a  very  liberal  fashion.  I  particularly 
appreciate  Jeffs  decision  to  hire  me  as  a  Draper  Fellow! 

Thanks  are  due  to  many  others.  Edward  McCormack  and  his  machine  shop  crew 
fabricated  parts  for  my  experiments.  Ken  Kaiser  consulted  (and  furnished  equipment)  for 
temperature  control  issues.  Dick  Caruso  quickly  processed  purchase  requests.  Mert  Prince, 
Connie  Cardoso,  and  Brenda  Hugh  helped  me  Hbricate  samples,  often  at  the  last  minute. 
George  Mitchell  and  Kevin  Hoxie  helped  with  making  my  pointer  beam  hot  plate.  Greg 
Kirkos  consulted  on  mechanics  issue.  John  Kaupinnen  and  Ralph  Shaver  helped  get  me  set 
up  in  a  lab  for  the  sink/float  work.  Lance  Niles  provided  me  device  schematics.  The 
irradiations,  dosimetry,  and  fluence  estimates  were  done  by  Gordon  Kohse  at  the  MIT 
reactor,  with  Mary  Montesalvo  and  David  McKee  supporting  the  irradiations  at  UMass 
Lowell.  I  appreciate  the  criticism  (mostly  constructive)  and  suggestions  of  the  "radiation 
meeting"  regulars  and  guests:  Kaplesh  Kumar,  Ralph  Hopkins,  Nicole  Gerrish,  Andrew 
Meulenberg,  and  William  Kelleher. 

My  wife,  Sophie,  has  accepted  constraints  on  her  career  choices  that  result  from  my 
participation  in  the  Draper  Fellow  program  and  service  in  the  Air  Force.  She  has  also 
supported  me  both  emotionally  and  logisticaUy  during  the  more  demanding  phases  of 
classes,  qualifying  exams,  and  research.  I  owe  a  very  great  deal  to  her  advice,  tolerance,  and 
optimism,  as  well  as  to  the  example  of  discipline  and  courage  she  set  for  me  while  she 
earned  her  own  graduate  degree. 


This  thesis  was  prepared  at  The  Charles  Stark  Draper  Laboratory,  Inc.,  under  IR&D  project 
number  1002. 


Publication  of  this  thesis  docs  not  constitute  approval  by  Draper  or  the  sponsoring  agency 
of  the  findings  or  conclusions  contained  herein.  It  is  published  for  the  exchange  and 
stimulation  of  ideas. 


Clark  L.  Allred 


4 


THE  VIEWS  EXPRESSED  IN  THIS  ARTICLE 
ARE  THOSE  OF  THE  AUTHOR  AND  DO  NOT 
REFLECT  THE  OFFICIAL  POLICY  OR 
POSITION  OF  THE  UNITED  STATES, 
DEPARTMENT  OF  DEFENSE,  OR  THE  U.S. 

GOVERNMENT 


ASSIGNMENT 


Draper  Laboratory  Report  Number  T-1383 


In  consideration  for  the  research  opportunity  and  permission  to  prepare  my  thesis  by  and  at 
The  Charles  Stark  Draper  Laboratory,  Inc.,  I  hereby  assign  my  copyright  of  the  thesis  to  The 
Charles  Stark  Draper  Laboratory,  Inc.,  Cambridge,  Massachusetts. 


Clark  L.  Allred  August  4,  2000 


5 


CONTENTS 

Page 


ABSTRACT  ^ 

ac:knowi.edgements  4 

CONTENTS  6 

LIST  Ol'  I'ICiURES  8 

IJSTOl' TABLES  10 

1  INTRODUCTION  H 

1.1  Scope  of  thesis  11 

1.2  Literature  Sun^ey  12 

1.3  Monte-Carlo  simulations  13 

1 .4  Gamma  irradiation  of  MEMS  pointer  beam  strain  gauges  14 

1.5  Neutron  irradiation  of  glass  samples  16 

2  BACKGROUND  10 

2.1  7\  physical  perspective  on  radiation  damage  in  glass  19 

2.2  Gamma  radiation  damage  24 

2.3  Neutron  radiation  damage  27 

2.4  Glass  compostions  30 

3  MONTE  CARLO  SIMULATIONS  33 

3.1  Introduction  to  TRIM  33 

3.2  Glass  models  35 

3.3  Fast  neutron  damage  simulations  37 

3.4  Thermal  neutron  damage  simulations  43 

4  GAMMA  IRRADIATION  OF  POINTER  BEAMS  47 

4.1  Pointer  beam  device  47 

4.2  Experimental  procedure  49 

4.3  Results  35 

4.4  Discussion  61 

5  NEUTRON  IRRADIATION  OF  GLASS  65 

5.1  Overview  65 

5.2  Experimental  Procedure  65 

5.3  Results  and  Discussion  ^6 

6  DISCUSSION  85 


6 


TABLE  OF  CONTENTS  (continued) 

7  CONCLUSION  90 

7.1  Summary  of  principal  findings  90 

7.2  Importance  of  findings  to  MEMS  93 

7.3  Suggestions  for  fumre  work  94 

APPENDIX  A:  Pointer  beam  strain  measurements  95 

APPENDIX  B:  Sink/ float  measurements  99 

REFERENCES  101 


7 


LIST  OF  FIGURES 

Page 

2.1  Silica  tetrahedron  and  borate  triangle  used  in  CRN  model  19 

2.2  Continuous  random  network  20 

2.3  Radiation  damage  processes  in  solids  22 

2.4  Pyrex  strain  vs.  absorbed  gamma  dose  25 

2.5  Pyrex  strain  vs.  absorbed  gamma  dose,  low  dose  regime  26 

2.6  Strain  in  Pyrex  and  siUca  vs.  absorbed  gamma  dose  28 

2.7  Thermal  neutron-incuded  strain  29 

2.8  Fast  neutron-induced  strain  in  vitreous  silica  31 

3.1  Schematic  of  1  MeV  neutron  collision  cascade  in  Pyrex  34 

3.2  TRIM-calculated  vacancies  in  oxygen  recoil  cascade  42 

3.3  Schematic  thermal  neutron  capture  and  resulting  collision  cascade  45 

4.1  Pointer  beam  strain  gauge  47 

4.2  Pointer  beam  array  on  glass  wafer  50 

4.3  Pointer  beam  holder  51 

4.4  Placement  of  dosimeters  52 

4.5  Dose  rate  at  the  sample  holder  53 

4.6  Gamma-induced  strain  in  a  Pyrex  glass  wafer  56 

4.7  Gamma-induced  strain  in  Hoya  SD-2  glass  wafers  57 

4.8  Gamma-induced  strain  in  Pyrex  glass  wafer,  alternate  ^'^to  orientation  59 

4.9  Gamma-induced  strain  in  Hoya  SD-2  glass  wafer,  alternate  ^"^Co  60 

orientation 

5.1  Cut  Pyrex  wafer  and  sample  assignments  66 

5.2  Glass  sample  pair  67 


8 


5.3  Sink/float  apparatus  schematic  71 

5.4  Neutron  irradiation-induced  strain  in  Pyrex  as  measured  by  sink/ float  77 

5.5  Neutron  irradiation-induced  strain  in  Hoya  SD-2  as  measured  by  78 

sink/  float 

5.6  Neutron  irradiation-induced  strain  in  Pyrex  as  measured  by  sink/ float  80 

(low-dose  regime) 

5.7  Neutron  irradiation-induced  strain  in  Hoya  SD-2  as  measured  by  81 

sink/ float  (low-dose  regime) 

5.8  Neutron  irradiation-induced  strain  in  Pyrex  and  Hoya  SD-2  as  82 

measured  by  sink/ float 

6.1  Comparison  of  tliermal  neutron-induced  strain  results  86 

6.2  Comparison  of  pointer  beam  to  previous  sink/ float  work  88 

7.1  Thermal  neutron-induced  strain  91 

7.2  Fast  neutron-induced  strain  92 


9 


T  JST  OF  TABLES 


Page 

1.1  TRIM  calculation  results  ^4 

1.2  Pointer  beam  gamma  irradiation  results  15 

1.3  Neutron  irradiation-induced  strains  16 

2.1  XPS  results  for  composition  of  Pyrex  and  Hoya  SD-2  glass  32 

3.1  Glass  compositions  used  in  TRIM  calculations  36 

3.2  Primary  knock-on  atoms  input  into  TRIM  40 

3.3  Fast  neutron  TRIM  results  43 

3.4  Thermal  neutron  TRIM  results  46 

4.1  Pointer  beam  pre-irradiation  positions  50 

4.2  Dosimeter  readings  54 

4.3  Pointer  beam  gamma  irradiation  results  58 

4.4  Pointer  beam  gamma  irradiation  results  (alternate  orientation)  58 

4.5  Pointer  beam  irradiation  averaged  results  61 

5.1  Glass  sample  irradiation  schedule  68 

5.2  Average  scatter  of  six  measurements  and  resultant  uncertainty  of  mean  76 

5.3  Linear  fitting  results  for  sink/ float  data  83 

5.4  Strain  in  glass  related  to  neutron  fluence  and  displacements  per  atom  83 

7.1  Neutron  irradiation-induced  strains  90 

7.2  Pointer  beam  gamma  irradiation  results  93 


10 


1  INTRODUCTION  AND  OVERVIEW 


1.1  Scope  of  thesis 

The  goal  of  the  research  reported  in  this  thesis  is  to  understand  dimensional  changes 
that  a  borosilicate  glass  substrate  can  undergo  when  exposed  to  radiation.  In  order  for  glass 
substrate-mounted  MEMS  devices  (such  as  silicon  oscillating  beams)  to  behave  predictably 
when  used  in  a  radiation  environment,  one  must  understand  the  effect  of  that  environment 
on  material  properties  such  as  density,  coefficient  of  thermal  expansion,  and  elasticity  for  the 
materials  involved  in  the  system.  This  study  treats  only  density  (and  hence,  length)  changes 
which  can  set  up  differential  strains  between  an  oscillator  and  its  substrate,  affecting  the 
stress-dependent  frequency  of  beam  oscillation. 

Forms  of  radiation  tliat  are  considered  here  include  1  MeV  neutrons  (referred  to 
below  simply  as  ’Tast”  neutrons),  thermal  neutrons,  and  1  MeV  gamma  rays.  Real  radiation 
environments  (^.^.,  outer  space,  nuclear  reactors)  wiQ,  of  course,  have  a  broader  energy 
spectrum  for  neutrons  and  gamma  rays.  Typical  nuclear  reactor  environments  include 
neutron  energy  spectra  with  centroid  energies  not  far  from  1  MeV.  For  the  average 
energy  of  emitted  neutrons  is  about  1.68  MeV.^  ^  Shielding  between  fissioning  nuclei  and  the 
target  can  modify  this  average  energy.  For  precise  radiation  damage  calculations,  geometry 
and  nature  of  surrounding  materials,  angle  of  incidence,  and  incident  energy  spectrum  for 
the  radiation  must  be  known.  It  is  nevertheless  useful  to  model  a  simpler  situation 
(monoenergetic  neutrons  and  gamma  rays)  to  obtain  order-of-magmtude  predictions  of 
radiation-induced  dimensional  changes.  These  predictions  can  then  be  scaled  by  taking  into 
account  the  amount  of  energy  (relative  to  1  MeV)  deposited  in  a  material  by  radiation 
particles. 
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Two  glasses  chosen  principally  for  their  close  thermal  expansion  coefficient  match 
with  silicon,  Pyrcx®'  “  and  Hoya  SD-2^^^^,  were  studied.  The  compaction  rates  of  Hoya  SD-2 
and  Pyrex  under  fast  neutron  bombardment  were  established  for  the  first  time.  The 
compaction  rate  of  Hoya  SD-2  with  thermal  neutron  fluence  was  measured  for  the  first  time, 
and  the  larger  compaction  rate  of  Pyrex  witli  thermal  neutron  fluence  was  confirmed. 
Evidence  for  increased  radiation-resistance  in  anodically  bonded  Pyrex  and/ or  radiation- 
induced  creep  at  the  anodic  bond  interface  was  found  with  the  novel  application  of  a  MEMS 
pointer  beam  strain  gauge  in  a  gamma  irradiation  experiment.  It  was  shown  that  the 
differential  strain  induced  by  gamma  radiation  in  MEMS  devices  anodically  bonded  to  glass 
wafers  is  much  less  in  magnitude  than  would  be  expected  by  considering  the  strains  that 
would  have  been  caused  by  similarly  irradiating  the  raw  glass  and  silicon. 

The  approach  taken  here  to  the  problem  can  be  summarized  in  four  steps:  1)  a 
survey  of  the  literature,  2)  Monte-Carlo  simulations  of  thermal  and  fast  neutron  collisions 
cascades,  3)  gamma  ray  irradiations  of  both  glasses  using  MEMS  strain  gauges  to  measure 
expansion/compaction,  and  4)  neutron  irradiation  of  both  glasses  followed  by  a  floatation 
technique  ("sink/float”)  used  to  measure  density  changes. 

1.2  Literature  survey 

Compaction  of  Pyrex  with  gamma  irradiation  has  been  studied  by  Shelby^  \ 

Conners^  \  Sato^  ^  and  Zdaniewski.^  ^  Of  these  studies,  Shelby’s  covers  the  greatest  dose 
range.  Since  Conners’  and  Sato’s  finding  support  Shelby’s  findings,  Shelby’s  work  is  used 
here  as  a  reference  for  our  findings.  Shelby's  data  reveal  an  initially  linear  trend  of 
compaction  with  dose  of  -1.4  X  10  strain  per  rad  at  doses  up  to  ^^2  X  10  rad.  The  effect 
begins  to  saturate  at  higher  doses. 
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The  effects  of  thermal  neutron  irradiation  in  Pyrex  were  studied  by  Paymal.^  "  He 
employed  a  floatation  technique  called  "sink/ float"  (the  same  technique  used  later  by  Shelby) 
to  measure  small  irradiation-induced  density  differences  in  several  glass  types.  He  observed 
a  compaction  rate  with  dose  that  at  lower  thermal  neutron  fluences  (<  1  X  10^^  n/crn)  was 
approximately  linear,  with  a  slope  of  -2  X  10  strain  per  n/crn. 

No  studies  involving  fast  neutron-induced  strain  are  known  for  either  glass. 

Primak/  however,  did  report  the  initial  compaction  rate  of  vitreous  silica  with  fast  neutron 
exposure  to  be  -5.6  X  10  strain  per  n/ cm^  Shelby’s^  results  for  gamma  irradiation  of 
vitreous  silica  show  that  Si02  compacts  an  order  of  magnitude  less  quickly  with  dose  than 
does  Pyrex. 

1.3  Monte-Carlo  simulations 

One  mechanism  for  tlie  causing  dimensional  changes  in  glass  is  the  collision  of  a  fast 
neutron  with  atoms  in  the  glass  network  structure.  The  atoms  so  struck  can  be  forced  from 
their  network  sites  with  substantial  kinetic  energy,  in  turn  leading  to  other  atomic 
displacements  and  subsequent  network  rearrangements.  This  sequence  of  atomic 
displacements  is  referred  to  as  a  collision  cascade,  and  the  atom  directly  displaced  by  an 
incident  fast  neutron  is  called  the  primary  knock-on  atom  (PICA). 

Collision  cascades  can  also  be  caused  by  nuclear  reactions  that  yield  energetic 
particles.  The  reaction  of  interest  in  this  thesis  is  the  absorption  of  a  thermal  neutron  by 
an  isotope  20%  naturally  abundant  in  boron  and  present  in  both  Pyrex  and  Hoya  SD-2. 
Upon  absorbing  a  thermal  neutron,  the  nucleus  splits  into  a  lithium  nucleus  with  a 
kinetic  energy  of  ~0.9  MeV,  and  an  alpha  particle  with  kinetic  energy  of  ~  1.6  MeV.  Each  of 
these  products  can  then  each  initiate  a  collision  cascade. 
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Both  thermal  and  fast  neutron  colUsion  cascades  involve  similar  PfCA  energies,  and 
the  majority  of  displaced  atoms  will  be  displaced  by  ions  created  from  secondary  coUisions. 
Therefore,  it  seems  reasonable  to  assume  that  the  network  rearrangement  processes  involved 
in  both  types  of  collision  cascades  will  be  similar,  and  that  the  dimensional  change  in  a  glass 
can  be  related  directly  to  the  number  of  atoms  displaced  in  these  collision  cascades.  To 
relate  fast  neutron-induced  compaction  to  thermal  neutron-induced  compaction,  we  need  an 
idea  of  how  many  displacements  occur  in  each  type  of  cascade. 

Monte-Carlo  simulations  using  TRIM  (version  2000.38),  a  computer  program  widely 
used  in  the  ion-implantation  industr)',  were  carried  out  for  both  thermal  and  fast  neutron 
collision  cascades  in  both  Pyrex  and  Hoya  SD-2.  Inputs  to  TRIM  included  glass 
composition  and  the  energy  distribution  of  PKAs.  Table  1.1  summarizes  the  results  for 
glass  wafers  780  pm  thick. 


Pyrex 

Hoya  SD-2 

Displacements  per  incident 
fast  neutron 

22 

24 

Displacements  per  incident 
thermal  neutron 

100 

31 

Table  1.1  TRIM  calculation  results 


1.4  Gamma  irradiation  of  MEMS  pointer  beam  strain  gauges 

A  MEMS  pointer  beam  strain  gauge,  developed  at  Charles  Stark  Draper  Laboratory 
(CSDL)  to  monitor  differential  strain  in  the  silicon  pointer  beam  and  the  glass  substrate  to 
which  it  is  anodically  bonded,  was  employed  in  a  gamma  irradiation  experiment.  The  gauges 
(referred  to  hereafter  as  "pointer  beams")  consist  of  a  beam  suspended  by  and  perpendicular 
to  two  beams  which  arc  slightly  offset  where  tiiey  join  the  main  beam.  These  two  beams  are 
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in  turn  suspended  from  silicon  columns  that  are  anchored  to  the  glass.  The  two  support 
beams  are  offset  so  that,  when  the  anchor  points  arc  moved  further  apart  or  closer  together, 
the  main  beam  is  caused  to  rotate  slightly.  Teeth  at  the  end  of  tlie  main  beams  serve  as 
reference  marks  for  beam  movement. 

The  gamma  radiation  was  expected  to  cause  dimensional  changes  in  the  glass,  but 
not  in  the  silicon.^  Any  pointer  beam  displacement,  it  was  assumed,  should  correspond 
direedy  to  strain  in  the  glass  substrate.  The  slopes  of  Hnear  fits  to  the  data  obtained  are 
shown  in  Table  1.2.  The  uncertainty  indicated  in  the  slope  is  due  to  two  alternative  gamma 
source  geometries  of  which  we  were  unaware  at  the  time  of  the  experiment.  Reduced 
values  for  the  fits  were  close  to  0.7  for  Pyrex  and  0.5  Hoya  SD-2,  irrespective  of  source 
geometry.  This  indicates  that  while  the  slope  on  the  fitted  lines  is  small,  it  is  of  statistical 
significance. 


Strain/ rad 

Pyrex 

5.24  (±  3.78)  X  10  ’^ 

Hoya  SD-2 

4.50  (±0.16)  X  10  '" 

Table  1.2  Pointer  beam  gamma  irradiation  results 

Shelby’s  results  for  Pyrex  provide  a  point  of  reference  for  these  results.  The  pointer 
beams  results  indicate  a  strain/rad  slope  26  times  smaller  than  Shelby  reports.  We  propose 
two  possible  explanations  for  this  discrepancy:  1)  the  anodic  bonding  process  has 
significandy  altered  the  radiation  sensitivity  of  the  Pyrex,  and/ or  2)  radiation-induced  creep 
occurred  at  the  anodic  bond  interface,  allowing  the  pointer  beam  anchor  points  to  slip  as  the 
glass  wafer  expanded.  These  results  are  of  great  sigmficance  to  the  operation  of  MEMS 
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devices,  since  they  indicate  a  smaller  net  radiation  response  for  anodically  bonded  MEMS 
devices  than  was  predicted  based  upon  the  responses  of  the  individual  materials. 

1.5  Neutron  irradiation  of  glass  samples 

Five  samples  of  each  glass  were  irradiated  with  neutrons  in  MIT's  reactor  and  then 
used  in  a  sink/ float  experiment  to  determine  irradiation-induced  density  changes.  Exposure 
times  ranged  from  30  s  to  6000  s  in  a  neutron  flux  comprised  of  1  X  10’ '  n/ eax  j s  fast 
neutrons  and  6.7  X  lO’^  n/cm^/s  thermal  neutrons.  Linear  fits  were  made  to  the  data  with 
reduced  values  of  1.01  for  Pyrex  and  1.58  for  Hoya  SD-2.  Compaction  was  apportioned 
between  thermal  and  fast  neutron  collision  cascades  using  the  displacement  numbers 
calculated  using  TRIM,  with  the  results  indicating  that  less  than  1%  of  the  compaction  was 
due  to  fast  neutrons.  Using  the  TRIM  displacement  numbers,  however,  one  may  make  a 
prediction  of  tire  response  of  the  glasses  to  both  tliermal  and  fast  neutrons.  The 
strain/ flucnce  relationships  thus  obtained  are  listed  in  Table  1.3,  along  with  values  relating 
strain  to  displacements  per  atom  (dpa). 


Strain/n/cm^ 

thermal 

Strain/n/cm^ 
fast  (1  MeV) 

Strain/dpa 

Pyrex 

-2.77xl0‘^" 

-6.07xl0‘^' 

-1.49 

Hoya  SD-2 

-1. 01x10^’ 

-7.89x10'^^ 

-0.19 

Table  1.3  Neutron  irradiation-induced  strains 


The  Hoya  SD-2  is  seen  to  undergo  almost  an  order  of  magnitude  less  compaction 
than  the  Pyrex  at  a  given  fast  neutron  fluence.  The  response  of  Hoya  SD-2  to  thermal 
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ncutxons  is  even  smaller  in  proportion  to  that  of  the  Pyrex;  this  can  be  partly  attributed  to 
the  smaller  boron  concentration  in  Hoya  SD-2. 

The  values  for  fast  neutron-induced  compaction  in  Pyrex  and  Hoya  SD-2  are  based 
upon  the  assumption  of  equivalent  damage  mechanisms  for  thermal  and  fast  neutron 
collision  cascades.  We  can  test  this  assumption  by  comparing  our  results  with  those  of  the 
studies  cited  in  section  1.2  and  check  for  consistency.  The  strain/ fluence  relationsliip  for 
Pyrex  is  in  good  agreement  with  Paymal's  result.  There  are  no  data  available  for  direct 
comparison  with  the  fast  neutron  results,  but  comparison  of  the  fast  neutron-induced  strain 
in  Pyrex  with  Primak’s  results  for  vitreous  silica  demonstrates  that  the  Pyrex  compacts  by 
about  an  order  of  magnitude  more  than  does  vitreous  silica  at  a  given  dose.  Shelby  showed 
that  this  same  relative  response  in  radiation  compaction  holds  for  gamma  irradiation, 
providing  (albeit  indirect)  support  for  validity  of  the  assumption  that  thermal  and  fast 
neutron  damage  can  be  compared  through  the  relative  number  of  displacements  in  their 
collision  cascades. 
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2  BACKGROUND 


2.1  A  physical  perspective  on  radiation  damage  in  glass 

2.1.1  Structure  of  glass 

Glass  is  a  metastable  solid  that  can  be  prepared  from  the  melt  by  rapid  cooling. 
"Rapid"  is  a  relative  term:  in  this  case,  it  means  rapid  enough  that  the  material  in  question 
does  not  have  enough  time  to  undergo  the  kinetic  processes  involved  in  crystallization.  For 
a  simple  crystal  structure,  such  as  a  close-packed  pure  metal,  this  is  very  rapid  indeed.  A 
metallic  glass  can  be  prepared  by  quenching  at  a  rate  of  about  10^’  °C/ s.^  ^  In  the  case  of  an 
oxide  such  as  silica,  the  atomic  organization  required  to  maintain  local  charge  neutrality 
makes  crj^stallization  a  slower  process,  and  vitreous  silica  can  be  prepared  from  the  melt 
simply  by  cooling  to  room  temperature  at  a  cooling  rate  of  about  °C/ s. 

A  model  for  the  structure  of  covalently  bonded  glasses  was  proposed  by 
Zachariasen^^.  Now  known  as  the  continuous-random-network  model,  it  describes  an  oxide 
glass  as  being  a  randomly  arranged  network  of  polyhedra  that  are  the  basic  oxide  units. 
Polyhedra  are  shown  in  Figure  2.1  for  Si02  and  B2O3,  two  key  components  of  borosilicate 
glass. 


B2O3 

Figure  2.1  Silica  tetrahedron  and  borate  triangle  used  in  random  continuous  network 
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These  polyhedra  can  be  joined  at  the  corners  (by  means  of  a  bridging  oxygen)  to 
form  a  network.  In  silica,  the  O-Si-O  bond  angles  are  rigidly  fixed  due  to  the  four- fold 
coordination  of  the  silicon  atoms,  but  the  oxygen  atoms  arc  only  bonded  to  two  silicon 
atoms.  A  greater  flexibility  in  the  A-O-A  bond  angles  results.  The  basic  polyhedra  can 
therefore  link  into  a  network  composed  of  rings  containing  varying  numbers  of  member 
polyhedra.  A  two-dimensional  representation  of  such  a  network  is  shown  in  Figure  2.2. 


Figure  2.2  Continuous  random  network 


The  addition  of  alkali  atoms  such  as  sodium  to  a  borosilicate  glass  alters  its  chemical 
structure  and  physical  properties.  Alkali  atoms  act  as  network  7nodijiers  in  silica,  reducing 
connectivity  in  the  network^^^  by  providing  local  charge  neutrality  for  non-bridging  oxygen 
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(an  oxygen  on  a  silica  tetrahedron  that  is  left  dangling).  When  and  are  present,  the 
addition  of  alkali  atoms  can  have  the  opposite  effect  of  increasing  connectivity  in  the 
network  by  providing  charge  neutrality  for  four-coordinated  A1  and  B  atoms  (B^^  and  Na^  ^ 
effectively  substituting  for  a  Si"^  ).  Physical  properties  that  change  with  the  addition  of  alkali 
atoms  include  increased  fracture  toughness^^  and  lowered  thermal  expansivity^^.  With 
increasing  alkali  content  (as  the  B/Na  or  Al/Na  ratio  falls  below  one),  connectivity  is  no 
longer  increased  and  Na^^  returns  to  its  role  as  a  network  modifier.  One  can  therefore 
expect  the  radiation  response  of  a  material  to  depend  not  only  on  A1  and  B  concentration, 
but  also  on  the  ratio  of  these  concentrations  to  tliat  of  alkali  atoms. 

2,1.2  Radiation-induced  alteration  of  glass  structure 

Radiation  has  several  mechanisms  for  affecting  a  material’s  structure.  Incident 
gamma  rays  can  interact  with  electrons  in  the  material  by  Compton  scattering.  Electrons  so 
scattered  can  then  baUistically  collide,  elastically  and  inelastically,  with  atoms  in  the  material. 
Both  gamma  rays  and  scattered  electrons  can  ionize  atoms,  as  well  as  give  rise  to  a 
phenomenon  called  radiolysis,  by  which  atomic  displacement  and  bond  rearrangement 
results  from  electronic  exitations^^.  Fast  neutrons  can  baUistically  collide  widi  atoms,  thus 
displacing  them  and  causing  related  electronic  defects.  Thermal  neutrons  can  provoke 
nuclear  reactions  (see  discussion  in  Chapter  3)  which  produce  energetic  ions  which  can  both 
electronicaUy  and  baUisticaUy  impart  energy  to  a  material.  A  summary  of  radiation 
interactions  in  materials  is  shown  in  Figure  2.3^^.  One  can  imagine  these  effects  in  the 
context  of  the  random  continuous  network  model.  Rupmred  bonds,  whether  caused  by 
radiolysis  or  collisional  removal  of  atoms  from  a  ring,  lead  to  atomic  rearrangements  and 
ring  size  changes.  These  ring-size  changes,  it  is  thought,  lead  to  density  changes. 
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THERIv'L\L  _ ABSORPTION^  EXCITED  _ decay  ^  IMPURITY 


igute  2.3  Radiation  damage  processes  in  solids 


Characterizations  in  the  literature  of  the  relationship  between  ring  size  (a  parameter 
that  radiation  may  change  in  a  material)  and  material  density  vary.  A  review  of  radiation 
effects  in  glasses^^'  cites  topological  modeling  studies^^  in  claiming  that  "larger  rings  are 
responsible  for  the  higher  density  of  the  denser  silica  polymorphs  because  they  can  fold  back 
on  themselves,  and  the  densification  of  vitreous  silica  by  irradiation  could  arise  from  opening 
up  of  6-membered  rings  to  larger  rings  in  regional  boundaries."  On  the  other  hand,  Ezz- 
Eldin^  ^  speculates  that  "damage  by  irradiating  species  can  cause  compaction  of  Si02  [silicate 
glasses]  by  breaking  bonds  between  tetrahedra,  allowing  the  formation  of  different  ring 
configurations  such  that  the  average  ring  size  is  smaller,  thus  leading  to  denser  structure." 

He  then  sites  evidence^  that  "smaller  rings  do  form."  A  more  recent  study^^^  states  that 
"whereas  it  might  be  supposed  that  the  presence  of  large  rings  implies  lower  density,  in  fact 
the  opposite  is  true...  As  the  polymoiphs  become  denser,...  the  ring  size  increases  and  the 
distribution  of  ring  size  widens."  The  author  goes  on  to  cite  tlie  example  of  crystalline  forms 
of  Si02,  and  indeed  finds  that  average  ring  size  trades  monotonically  with  density. 

Not  all  silicate  glasses  compact  under  irradiation.  Paymaf  reports  diat  while 
Pyrex  compacts  with  thermal  neutron  exposure,  a  borosiheate  glass  with  PbO  additions 
clearly  expands.  The  ring  model  says  nothing  about  whether  or  not  a  glass  will  expand  or 
compact  upon  radiation  —  it  only  relates  ring  size  to  density.  Shelby^^^'  wrote  that  "it 
becomes  quite  obvious  that  the  densification  mechanism  will  require  considerable  study 
before  the  details  of  the  process  become  fully  understood...  it  is  quite  possible  that  the 
results  reported  here  result  from  competing  processes  which  involve  both  positive  and 
negative  dilatations."  Paymaf  proposed  just  such  a  theor)^  for  his  neutron  damage 
obsen^ations  in  which  the  energy  deposited  by  fast  lithium  and  helium  ions  in  a  glass  is 
mostly  concentrated  in  a  small  end-of-trajectory  region,  causing  a  thermal  spike.  The  high- 
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temperature  melt  is  rapidly  cooled  by  the  glass  matrix,  "freezing  in"  an  expanded  region. 

The  glass  immediately  around  this  expanded  region  is  forced  to  compact.  The  relative 
volume  of  these  expanded/compacted  regions  determines  the  macroscopic  density  change. 
While  being  a  nice  conceptual  model,  this  theor)^  still  docs  not  address  the  fundamental  issue 
of  directly  relating  composition  to  radiation  response. 

There  is  evidence  that  boron  concentration  as  well  as  sodium  concentration  (relative 
to  that  of  boron)  affect  radiation  response.  Shelby^  writes  that  "it  is...  quite  possible  that 
the  effects  reported  here  arc  strongly  related  to  the  two-phase  nature  of  most  borosiheate 
glasses...  the  glass  consists  of  a  continuous  phase  with  a  composition  very  near  that  of 
vitreous  silica...  the  second  phase...  contains  most  of  the  alkali  and  boron  atoms.  This  alkali 
borate  phase  is  suspected  of  being  more  susceptible  to  radiation  compaction  than  vitreous 
silica,  with  the  overall  material's  compaction  depending  on  the  relative  volume  of  this  second 
phase,  as  well  as  the  coordination  of  the  boron  within  that  phase.  As  discussed  above,  this 
boron  coordination  will  depend  on  the  alkali  content  of  the  phase. 

2.2  Gamma  radiation  damage 

Shelby^  Conners^  Sato^  ^^  and  Zdaniewsky^  ^'^  report  measurements  of  gamma 
radiation-induced  density  changes  in  Pyrex.  All  four  employed  a  variant  of  the  sink/float 
method  described  in  Chapter  5.  Their  results  are  summarized  in  Figures  2.4  and  2.5. 

Ciood  agreement  is  found  between  results  of  Shelby,  Conners,  and  Sato.  The  mean 
difference  of  29  %  in  strain  between  the  results  of  Shelby  and  Conners  at  high  doses  may  be 
due  to  differences  in  glass  compositions  and/or  processing  history  that  may  have  existed 
between  glass  samples.  Zdaniewsky's  results  are  for  doses  of  10 , 10 , 10  ,  and  10  rad.  The 
compaction  he  measured  occurs  much  more  quickly  with  dose  than  in  the  other  studies. 
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Unfortunately,  while  citing  Shelby’s  study,  he  docs  not  comment  on  the  possible  reasons  for 
the  discrepancy  in  densification  trend. 

Shelby's  study  compared  the  compaction  of  several  borosilicate  glasses  to  that  of 
vitreous  silica.  He  concluded  that  with  increased  boron  content,  increased  compaction 
occurred.  A  comparison  of  strain  in  Pyrex  to  that  in  Si02  is  shown  in  Figure  2.6.  While 
having  the  same  slope  on  this  log-log  plot,  the  curve  for  Pyrex  shows  an  order  of  magnitude 
more  strain  at  a  given  gamma  dose.  It  should  be  noted  that  the  slope  of  compaction  vs. 
gamma  dose  (on  a  log-log  plot)  for  silica  varies  by  as  much  as  a  factor  of  three  in  studies  by 
Shelby^  Higby^^’’,  and  Primak^^\  Higby  attributes  this  variation  to  minor  differences  in 
impurities  or  in  the  degree  of  phase  separation.  It  is  not  surprising,  given  these  differences 
in  silica  results,  to  see  the  relatively  minor  variations  in  the  behavior  of  different  samples  of 
Pyrex  shown  in  Figure  2.4. 

23  Neutron  radiation  damage 

Paymal^  ^^  ^^^  carried  out  an  extensive  study  on  tire  effects  of  thermal  neutron  damage 
(the  mechanisms  for  which  are  discussed  in  Chapter  3)  on  the  properties  of  borosilicate 
glass.  His  model  of  competing  and  interdependent  regions  of  compaction  and  swelling, 
mentioned  above,  was  fitted  to  his  data  for  compaction  in  Pyrex.  The  resulting  relationship 
is: 

(2.1)  Ap  =  ^ 

In  Equation  2.1,  Ap  is  the  change  in  density  in  units  of  lO  '*  g/ cm’  and  D  is  the  thermal 
neutron  fluence  in  n/ctn.  A  conversion  of  this  relationship  for  strain  in  Pyrex  is  shown  in 
Figure  2.7. 
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Gamma  Dose,  rad 

absorbed  gamma  dose  (from  Shelby 


Paymal  attributes  the  saturation  of  strain  evident  in  Figure  2.7  to  overlap  of  the  end- 
of-trajectory  thermal  spike  regions.  A  similar  saturation  is  displayed  by  fast  neutron  damage 
in  vitreous  silica  in  results  from  Primak^^“.  These  results  are  shown  in  Figure  2.8.  The  shape 
of  this  cuiwc  agrees  with  that  of  the  Paymal  thermal  neutron  curve  for  Pyrex,  but  does  not 
reveal  the  relative  sensitivity  of  the  two  materials  to  atomic  displacements.  While  vitreous 
silica  will  be  unaffected  by  thermal  neutrons  (it  contains  no  boron),  it  is  nevertheless 
possible  to  compare  the  thermal  neutron  damage  in  Pyrex  to  the  fast  neutron  damage  in 
vitreous  silica  by  means  of  Monte  Carlo  simulations  (discussed  in  detail  in  Chapter  3).  These 
simulations  predict  22  atomic  displacements  in  Pyrex  per  incident  fast  neutron  and  100 
displacements  for  each  incident  thermal  neutron.  Results  will  be  slightly  different  for 
vitreous  silica,  but  with  similar  densities  for  Pyrex  and  silica  the  difference  will  be  minimal. 

If  we  consider  the  low-dose  linear  region  in  Figure  2.7  and  2.8  at  the  strain  of -1  X  10  ,  and 
multiply  the  corresponding  fluences  by  the  above  displacement  numbers  for  incident  fast 
and  thermal  neutrons,  we  obtain  that  to  achieve  this  strain  ^4  X  10^^  atomic  displacements 
arc  required  in  Pyrex,  while  -  4  X  10^'^  atomic  displacements  are  required  in  vitreous  silica. 
This  order-of-magnitude  difference  in  radiation  response  is  the  same  as  observed  for  gamma 
radiation  in  Figure  2.6. 

2.4  Glass  compositions 

Table  2.1  shows  compositions  for  the  glasses  under  consideration  here.  The  Evans 
I7ast^^^  company  was  employed  to  perform  X-ray  photoelectron  spectroscopy  (XPS)  analysis 
on  Pyrex®^^'*,  Borofloat®’ ^  (a  Pyrex  equivalent),  and  Hoya  SD-2®^^*  glass  wafer  samples. 
The  Pyrex  and  Borofloat  samples  were  found  to  be  compositionally  similar  to  within  the 
signal-to-noise  limitations  of  the  measurement,  and  so  have  been  averaged  together.  The 
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Fast  Neutrons  Fluence,  n/cm 

Figure  2.8  Fast  neutron-induced  strain  in  vitreous  silica  (from  Primak^^^) 


O 

Si 

B 

A1 

Na 

K 

Ca 

c 

Mg 

Zn 

Pyrex 
(this  study) 

63.6 

26.3 

6.4 

1.1 

1.2 

0.3 

- 

1.2 

- 

- 

Pyrex, 

Shelby^'' 

64.1 

26.9 

4.8 

0.8 

2.8 

0.4 

0.1 

- 

0.1 

“ 

Pyrex, 

PaymaP^ 

64.3 

27.6 

4.6 

0.5 

2.2 

0.7 

- 

-- 

- 

Hoya  SD-2 

61.0 

21.1 

1.1 

12.2 

0.5 

- 

- 

0.6 

3.0 

0.8 

Tabic  2.1  XPS  results  for  composition  (atomic  %)  of  Pyrex  and  Hoya  SD-2  glasses 


Hoya  SD-2  data  arc  also  the  average  of  two  measurements.  The  manufacturer-provided  data 
for  Pyrex  and  Hoya  SD-2  hst  densities  of  2.23  g/ cm^  and  2.60  g/ cm^,  respectively. 

It  can  be  seen  that  while  Pyrex  is  clearly  a  borosilicate  glass,  Hoya  SD-2  is  primarily 
an  aluminosilicate  glass,  modified  by  boron  and  magnesium.  Atom  number  densities  (based 
on  composition  and  mass  density)  are  found  to  be  approximately  6.90  X  10^^  / cm^  for  Pyrex 
and  7.65  X  10“  /cm’  for  Hoya  SD-2. 
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3  MONTE-CARLO  SIMULATIONS 


3.1  Introduction  to  TRIM 

Energetic  particles  incident  on  a  material  can  lose  energy  through  collisions  with 
atomic  nuclei  and  by  interactions  with  electrons.  Collisions  with  nuclei  can  result  in  atoms 
being  displaced  from  their  initial  positions  and,  then  as  energetic  ions,  colliding  with  other 
nuclei  and  causing  electronic  excitations  in  the  material.  This  sequence  of  primary, 
secondary,  and  higher-order  collisions  is  called  a  '^collision  cascade".  Figure  3.1  shows  a 
schematic  for  a  collision  cascade  resulting  from  a  1  MeV  neutron  striking  an  oxygen  atom  in 
Pyrex.  The  primar}^  knock-on  oxygen  atom  strikes  another  oxygen,  dislodging  it  from  its 
site,  then  both  atoms  continue  on  through  the  material,  diplacing  more  atoms  of  various 
species. 

The  Transport  of  Ions  in  Matter,  or  TRIM,  is  a  computer  code  available  free  of 
charge  from  James  F.  Ziegler  at  IBM  Research  as  part  of  the  SRIM  (Stopping  Range  of  Ions 
in  Matter)  software  package^  \  SRIM  runs  on  a  PC  with  a  DOS  operating  system.  Widely 
used  in  the  field  of  ion  implantation,  TRIM  performs  Monte-Carlo  simulations  of  collision 
cascades  for  any  ion  passing  through  elemental  materials  as  well  as  compounds.  Data 
provided  by  TRIM  includes  the  average  number  of  atomic  displacements  in  a  cascade,  the 
percentage  of  energy  lost  by  both  incident  ions  and  recoil  atoms  to  ionization  and  atomic 
displacement,  and  the  range  of  incident  ions  in  a  material.  One  can  even  generate  a  file 
containing  the  locations  of  all  collisions  within  a  cascade,  along  with  the  corresponding 
recoil  energies. 
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Schematic  of  1  MeV  neutron  collision  cascade 


TRIM  (vetvsion  2000.38)  simulations  were  carried  out  for  the  interaction  of  fast  and 
thermal  neutrons  in  Pyrex  and  Hoya  SD-2  glasses.  The  purpose  of  doing  this  was  two-fold: 

1)  data  were  found  in  the  literature'^  ^  for  thermal  neutron  damage  in  Pyrex,  and  an 
equivalence  was  needed  between  fast  and  thermal  neutron  collision  cascades  in  order  to 
estimate  fast  neutron  damage,  and  2)  planned  experiments  on  glass  would  involve  exposure 
to  both  thermal  and  fast  neutrons,  and  a  way  was  needed  to  apportion  the  observed  damage 
between  these  two  causes. 

The  key  piece  of  data  that  TRIM  can  furnish  to  help  estimate  this  fast 
neutron/ thermal  neutron  equivalence  is  the  number  of  displacements  per  collision  cascade. 

In  the  TRIM  lexicon,  "vacancies"  are  equivalent  to  the  "displacements"  referred  to  here.  In 
calculating  vacancies,  TRIM  takes  into  account  that  an  energetic  ion  that  transfers  enough 
energy  to  a  stationar}^  ion  to  dislodge  it  from  its  site  may  not  have  enough  energy  itself  after 
the  collision  to  escape  the  site.  In  this  case,  the  incoming  ion  replaces  the  ion  initially  at  the 
site,  and  no  vacancy  results.  "Displacements"  calculated  here  are  this  net  vacancy  result 
from  TRIM,  with  replacement  collisions  already  taken  into  account. 

3.2  Glass  models 

The  stopping  power  of  a  material  is  the  amount  of  energy  lost  to  that  material  by  an 
effluent  charged  particle  per  distance  travelled^^'.  Quoting  Ziegler,^^  "...the  stopping  [power] 
of  a  compound  may  be  estimated  by  the  linear  combination  of  the  stopping  powers  of  the 
individual  elements.  This  [Bragg’s]  rule  is  reasonably  accurate,  and  measured  stopping  of 
ions  in  compounds  usually  deviates  less  than  20%  from  that  predicted  by  Bragg's  rule." 
Furthermore,  "the  accuracy  of  Bragg's  rule  is  Umited  because  the  energy  loss  to  the  electrons 
in  any  material  depends  on  the  detailed  orbital  and  excitation  structure  of  the  matter,  and  any 
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differences  between  elemental  materials  and  compounds  will  cause  Bragg's  rule  to  become 


inaccurate." 

TRIM  accommodates  the  siinulation  of  ion  transport  in  compounds  like  glass  using 
a  "cores  and  bonds"  (OYB)  model.^^'  The  "cores"  are  closed-shell  atoms,  and  tlie  "bonds' 
are  the  bonding  valence  electrons  which,  depending  upon  the  exact  nature  of  a  bond,  will 
contribute  differing  amounts  to  the  overall  stopping  power  of  the  material  The  TRIM  file 
COMPOUND.DAT  contains  compositions  for  many  common  nuclear  materials,  and 
detailed  bond  information  on  a  subset  of  these. 


Element 

Pyrex 
atom  % 

Hoya  SD-2 
atom  % 

o 

65.0 

60.95 

Si 

25.0 

21.1 

B 

7.0 

1.05 

A1 

1.0 

12.2 

Na 

2.0 

0.5 

c 

- 

0.6 

Zn 

- 

0.75 

Mg 

- 

2.95 

Table  3.1  Glass  compositions  used  in  TRIM  calculations 


Table  3.1  shows  the  glass  compositions  used  for  the  TRIM  simulations.  The  Pyrex 
composition  differs  slightly  from  that  actually  measured  by  XPS  (see  Chapter  2).  There  are 
two  reasons  for  tliis:  1)  at  the  time  the  Pyrex  TRIM  calculations  were  performed,  the  XPS 
data  were  not  available  and  2)  a  model  for  Pyrex  was  already  included  in  the  TRIM  database. 
It  was  assumed  at  the  time  of  the  calculations  that  the  fuU  CAB  model  was  used,  but  further 
inspection  of  the  TRIM  COMPOUND.DAT  file  revealed  that  no  bonding  information  is 
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present.  It  is  therefore  presumed  that  Bragg's  rule  alone  is  being  used  in  the  TRIM 
calculations  for  both  Pyrex  (using  the  default  Pyrex  model)  and  Hoya  (with  manually  entered 
compositions). 

3.3  Fast  neutron  damage  simulations 
3.3.1  Introduction  to  fast  neutron  scattering 

TRIM  calculations  were  carried  out  for  1  MeV  neutrons  incident  on  both  Pyrex  and 
Hoya  SD-2  glass.  As  TRIM  will  not  accept  a  neutron  as  an  incident  particle,  the  nuclei 
recoiling  from  a  fast  neutron  collision  must  serv^e  this  purpose.  The  only  such  particle 
considered  here  was  the  first  atom  struck  by  a  fast  neutron,  or  the  primary  knock-on  atom 
(PKA).  For  energies  up  to  1  MeV,  the  cross  section  for  kinetic  energy  transferred  to  the 
nucleus  by  a  neutron  of  energy  can  be  approximated  by  ^  ^ : 

(3.1)  =  const/T,„ 
where 

(3.2)  0<T<T„  =  [4m, 1112/ (m,  +  m2)^]  E,  =  [4A/ (1+A)^]  E, 

In  the  above  equations,  m^  is  the  neutron  mass,  m2  is  the  target  atom  mass,  A  is  the  target 
atomic  mass  number,  and  T,^  is  the  maximum  transferred  kinetic  energy  possible.  Equation 

(3.2)  follows  from  classical  mechanics  using  conservation  of  energy  and  momentum  for  an 
incident  neutron  striking  a  nucleus  initially  at  rest.  Equation  (3.1)  says  that  all  possible  values 
for  transferred  kinetic  energy  are  equally  likely.  This  approximation  breaks  down  at  neutron 
energies  substantially  higher  than  1  MeV. 
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The  cross  section  for  scattering  of  a  fast  neutron  depends  upon  its  energy  as  well  as 
the  scattering  nucleus.  Neutron  scattering  cross  sections  are  available  in  the  literature  for  the 
naturally  occurring  isotopic  averages  of  elements  in  our  glasses/  ^  as  well  as  for  individual 
isotopes.  The  symbols  used  in  [3.9]  are  the  following: 

:  total  cross  section 

(T^.,  :  elastic  scattering  cross  section 

:  inelastic  scattering  cross  section 

^abs  •  absorption  cross  section 
a  :  nonelastic  cross  section, 

where:''’ 


^tot  “  ^el  ^inl  ^abs 


C^non  -  +  ^^abs 


For  1  MeV  incident  neutrons,  most  of  the  elements  in  Pyrex  and  Hoya  SD-2  have 
scattering  cross-sections  that  are  purely  elastic.  Exceptions  are  sodium  and  zinc,  which  have 
nonelastic  cross  sections  that  are  13%  and  4%,  respectively  of  their  elastic  cross  sections. 
These  small  nonelastic  contributions  to  scattering  are  made  even  smaller  for  the  overall  glass 
by  the  minor  presences  of  sodium  and  zinc  in  our  glasses.  Therefore,  the  simplify^ing 
assumption  of  elastic  scattering  was  made  for  all  primary  knock-ons. 

The  product  aN,  where  a  is  the  microscopic  cross  section  discussed  above  and  N  is 
die  number  density  of  atoms,  is  called  the  macroscopic  cross  section,  Z.  ^  The  mean  free 
path,  X,  of  a  neutron  is  given  by  1/S.  The  probability  of  a  neutron  not  scattering  before 
travelling  a  distance  x  in  a  material  is  given  by:'"’ 


(3.3) 


Pnot(x)  - 


-Ix 
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The  probability  of  a  neutron  scattering  before  travelling  a  distance  x  is  therefore: 

(3.4) 

3.3.2  TRIM  calculation  set-up 

TRIM  can  read  an  input  file  of  individual  recoiling  nuclei  and  calculate  a  collision 
cascade  for  each.  Since  our  glasses  contain  several  different  elements,  one  must  calculate  the 
relative  number  of  primary  knock-on  atoms  for  each  constituent  atom  type.  One  must  also 
account  for  different  energies  between  PI<j\’s,  as  given  by  Eq.  3.1. 

The  number  density  of  each  element  in  both  glasses  was  calculated  based  on  the 
density  of  the  glass  and  the  atomic  fraction  of  the  element.  Multiplying  tliis  number  density 
by  the  microscopic  cross  section  for  elastic  neutron  scattering  from  [3.6]  gives  the 
macroscopic  cross  section.  This  then  can  be  used  in  Eq.  3.4  to  give  the  fraction  of  incident 
neutrons  that  will  be  scattered  by  this  element.  The  glass  wafer  thickness  (the  distance  "x"  in 
Eq.  3.4)  used  for  these  calculations  was  780  pm,  that  of  a  typical  MEMS  substrate. 

The  most  probable  scatterer  in  Pyrex  and  Hoya  SD-2  is  oxygen,  scattering  2.8%  and 
2.9%,  respectively,  of  incident  1  MeV  neutrons.  Silicon  is  the  next  largest  contributor, 
scattering  around  0.5%  of  incident  neutrons  in  both  glasses,  with  other  elements  scattering 
relatively  smaller  numbers  of  neutrons.  Multiple  scattering  of  neutrons  is  therefore  ignored 
when  preparing  a  list  of  primary  knock-on  atoms  for  the  TRIM  calculation.  Each  scattering 
event  as  predicted  by  Eq.  3.4  is  considered  to  be  one  involving  a  1  MeV  neutron,  and  the 
total  probability  for  scattering  in  a  glass  is  assumed  to  be  the  sum  of  the  individual  scattering 
probabilities  of  each  element.  Furthermore,  the  neutron  flux  and  energy  is  assumed  to  be 
uniform  throughout  the  glass  wafer;  Le.^  no  consideration  is  given  to  the  fact  that  the  amount 
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of  energy  deposited  by  an  incident  neutron  flux  will  be  slightly  greater  near  the  entrance  side 
of  the  wafer  than  near  the  exit  side  due  to  attenuation  of  neutron  kinetic  energy.  The 
numbers  of  each  type  of  PKA  used  in  TRIM  calculations  arc  shown  in  Table  3.2. 


Element 

Pyrex 

Number  of  PKAs 

Hoya  SD-2 
Number  of  PKAs 

Energy  Range 
(keV) 

o 

1100 

1100 

( 0  -  221 ) 

Si 

262 

212 

(0-133) 

B 

23 

6 

(0-310) 

A1 

3 

65 

(0-138) 

Na 

14 

3 

(0-160) 

C 

3 

( 0  -  284 ) 

Zn 

__ 

6 

(  0  -  59  ) 

Mg 

— 

18 

(0-152) 

K 

3 

- 

(  0  -  97  ) 

Tabic  3.2  Primar)^  Knock-on  7\toms  input  into  TRIM 

We  note  that  the  composition  of  the  input  PICAs  varies  somewhat  from  the  TRIM 
default  Pyrex  composition  the  presence  of  potassium).  The  composition  used  to 
generate  the  values  in  Table  3.2  comes  from  PaymaP  \  Again,  the  TRIM  runs  for  Pyrex 
predate  the  XPS  measurements,  and  it  was  desirable  to  leave  the  native  TRIM  Pyrex  model 
intact  for  the  target  material.  The  schematic  shown  earlier  in  Figure  3.1  for  a  1  MeV  neutron 
collision  cascade  is  representative  of  Pyrex,  with  oxygen,  silicon,  and  boron  present  in 
correct  proportions  and  their  relative  si/^cs  indicating  their  cross  section  for  scattering  the 
fast  neutron  (ie.,  the  relative  probability  of  being  the  primary  knock-on  atom). 

The  motivation  for  running  proportional  numbers  of  PIXAs  was  to  reduce 
computation  time.  One  could  have  simply  run  thousands  of  each  type  of  PICA,  then 
averaged  the  results  together  proportionately  to  the  neutron  scattering  probability  of  each 
element,  as  described  above.  Instead  of  doing  this,  it  was  judged  that  ^1400  total  PKAs 
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gave  acceptable  statistical  certainty  in  the  averaged  result  for  displacements  per  PICA,  and 
this  total  number  was  shared  among  the  elements  according  to  their  importance.  The  1100 
oxygen  PKAs  took  about  6  hours  to  calculate  on  a  Pentium  120  MHz  computer. 

The  last  column  in  Table  3.2  shows  the  range  for  possible  energy  transfer  from  the 
incident  1  MeV  neutron  to  the  PIC\  (see  Eq.  3.2).  To  simulate  the  spectrum  of  possible 
starting  energies  of  primary  knock-on  atoms  (Eq.  3.1),  equally  spaced  energies  for  PKAs 
were  included  in  the  TRIM  input  file.  For  example,  100  of  the  1100  oxygen  PICAs  started 
with  0  keV,  100  others  started  with  22.1  keV  (being  one  tenth  of  T^^  —  221  keV),  and  so  on 
up  to  the  last  100  atoms  which  started  with  221  keV.  The  PKA  energies  for  the  other 
elements  were  similarly  represented  by  0,  T^/10,  ...,  T^.  For  elements  having  n  fewer  than 
eleven  PKAs,  the  energy  range  was  simply  broken  into  n  intervals.  TRIM  calculations 
carried  out  later  revealed  that  assigning  aU  PKyVs  an  energy  of  T^^.^  =  Tj^/2  yields  very  similar 
results  to  the  above  procedure,  />.,  the  number  of  displacements  in  a  collision  cascade  varies 
fairly  linearly  with  initial  PKA  energy  in  these  energy  ranges.  Figure  3.2  shows  the  locations 
in  the  X-Y  plane  of  displacements  calculated  by  TRIM  for  two  oxygen  recoil  coUision 
cascades.  One  is  for  the  case  of  maximum  energy  transfer  of  221  keV,  the  other  for  one- 
tenth  of  that. 
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Figure  3.2  TRIM-calculated  vacancies  in  oxygen  recoil  cascade  in  Pyrex 


The  results  of  the  fast  neutron  TRIM  calculations  are  summarized  in  Table  3.3. 


Pyrex 

Hoya  SD-2 

Displacements  per  fast 
neutron  collision 

614 

644 

Probability  of  fast  neutron 
collision 

0.036 

0.038 

Displacements  per  incident 
fast  neutron 

22 

24 

Table  3.3  Fast  neutron  (1  MeV)  TRIM  results 

3.4  Thermal  neutron  damage  simulations 

3.4.1  Introduction  to  boron  thermal  neutron  capture 

Incident  thermal  neutrons  do  not  pose  any  threat  of  coUisionally  displacing  atoms  in 
the  two  glasses  considered.  Their  0.025  eV  average  energy  is  well  below  the  15-30  eV 
displacement  threshold  for  atoms  in  the  glass  network.  We  are  instead  concerned  with  the 
recoiling  alpha  particle  (helium  ion)  and  the  recoiling  lithium  ion  that  result  from  the  capture 
of  a  thermal  neutron  by  a  nucleus.  The  reaction  is:^  ^ 

(3.5)  'sB+oH  ->  2^6+ 3Li  +  2.5MeV 

The  2.5  MeV  is  apportioned  between  the  He  and  Li  ions  according  to  their  masses: 
E„,  =  (2.5  MeV)  M^  /  (M„,,  +  MJ  =  1.59  MeV 
Ei_,  =  (2.5  MeV)  M„,  /  (M„,.  +  M,  ,)  =  0.91  MeV 


43 


TRIM  calculations  arc  therefore  simpler  to  set  up  for  the  thermal  neutron  case  than 
for  fast  neutrons,  as  one  has  simply  to  calculate  displacement  of  monoenergetic  helium  and 
lithium  ions,  then  add  the  two  results.  Figure  3.3  shows  a  schematic  for  a  thermal  neutron 
capture  by  a  boron  atom  in  Pyrex.  A  boron  atom,  upon  absorbing  the  incident  thermal 
neutron,  is  replaced  by  energetic  lithium  and  helium  particles. 

3.4.2  Thermal  neutron  TRIM  set-up  and  results 

The  cases  of  incident  helium  and  hthium  particles  were  calculated  separately  with 
TRIM.  The  sum  of  the  average  vacancies  per  collision  cascade  gives  the  effect  of  a  boron 
neutron  absorption.  Because  manually  created  files  can  allow  TRIM  to  simulate  ions 
originating  within  a  material  (as  would  be  the  case  for  neutron  absorption),  both  helium  and 
lithium  originations  were  considered  for  this  case.  The  results  were  the  same  for  ineident 
ions  and  for  ions  starting  from  inside  the  glass  (J.e.,  the  surface  binding  energy  information 
included  in  the  TRIM  calculation  had  ncghgible  effects  on  the  results).  Overnight  runs  were 
carried  out,  accumulating  ~  17,000  lithium  ions  and  ~38,000  helium  ions  (the  statistics  are 
stable  long  before  this  many  ions  are  treated). 

As  with  the  fast  neutron  case,  one  must  determine  how  many  of  the  incident  thermal 
neutrons  will  be  absorbed  in  the  glass.  The  cross  section  for  thermal  neutron  absorption  for 
’"B  is  3837  barns  (1  barn  =  1  X  10  ^'*  cm^),  and  that  of  namrally  occurring  boron  is  767 
barns.’ "  Using  this  cross  section,  the  number  density  of  boron  in  Pyrex  (3.2  X  10^’ 
boron/cm’)  and  Hoya  SD-2  (8.0  X  10^''  boron/cm’),  and  a  glass  thickness  of  780  pm  in  Eq. 
3.4  yields  probabilities  for  thermal  neutron  of  absorpdon  of  17.3%  and  4.7%  in  Pyrex  and 
Hoya  SD-2,  respectively. 

The  results  of  the  TRIM  calculations  are  summarized  in  Table  3.4. 
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Thermal  neutron 


Figure  3.3  Schematic  thermal  neutron  capture  and  resulting  collision  cascade 


Pyrex 

Hoya  SD-2 

Displacements  per  thermal 
neutron  absorption 

578 

664 

Probability  of  thermal 
neutron  absorption 

0.173 

0.047 

Displacements  per  incident 
thermal  neutron 

100 

31 

Tabic  3.4  Thermal  neutron  TRIM  results 


To  calculate  displacements  per  atom  for  the  values  in  Table  3.3  and  Table  3.4,  one 
may  use  the  atom  number  densities  (based  on  composition  and  mass  density)  of  6.90  X  10^^ 
/cm  for  Pyrex  and  7.65  X  10“^  /cm/  for  Hoya  SD-2. 
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4  GAMMA  IRRADIATION  OF  POINTER  BEAMS 
4.1  Pointer  beam  device 

A  MEMS  pointer  beam  strain  gauge  has  been  developed  at  Charles  Stark  Draper 
Laboratory  (CSDL).  Its  primary  use  has  been  as  a  monitoring  tool  for  wafer  bonding  and 
fabrication  processes,  as  it  gives  an  indication  of  differential  strains  present  between  a  glass 
substrate  and  the  silicon  wafer  bonded  to  it. 

A  diagram  of  tlie  pointer  beam  device  is  shown  in  figure  4.1. 


1  mm 


Figure  4.1  Pointer  beam  strain  gauge 

The  areas  shaded  in  black  represent  portions  of  the  device  that  are  anodically  bonded 
to  the  glass  substrate.  The  hashed  regions  represent  silicon  parts  suspended  above  the 
surface  of  the  substrate.  The  entire  suspended  portion  is  thus  supported  by  the  two  silicon 
blocks  at  the  bottom  of  Figure  4.1. 
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The  key  to  the  strain  gauge  lies  in  the  slight  offset  between  the  two  horizontal 
supports  for  the  vertical  beam.  Thus,  if  the  silicon  is  in  compression,  the  top  of  the  vertical 
beam  will  swing  to  the  left.  If  the  silicon  is  in  tension,  it  will  swing  to  the  right.  The  row  of 
teeth  spaced  at  5  pm  attached  to  the  vertical  beam,  when  referenced  to  tbe  row  of  teeth 
spaced  at  4.5  pm  on  the  bonded  silicon  block  at  the  top  of  Figure  4.1  serve  as  a  vernier 
caliper,  allowing  one  to  measure  the  degree  of  compression  or  tension.  The  beam  at  the 
bottom  of  Figure  4.1  that  runs  parallel  to  the  offset  beams  plays  no  role  in  the  strain 
measurement:  its  detached  end  (on  die  right)  seiwes  to  measure  strain  gradients  in  the  silicon 
as  it  curls  out  of  plane. 

A  finite  element  analysis  performed  by  CSDT  staff  showed  that  the  displacement  of 
the  pointer  beam  teeth  with  respect  to  the  anchored  reference  teeth  is  given  by: 

(4.1)  d  =  (0.035pm)£ , 

where  d  is  the  tip  displacement  in  microns  and  £  is  the  strain  in  parts  per  milhon  (ppm). 

Since  the  vernier  reads  to  1  /1 0  of  a  tooth  width,  or  0.5  pm,  this  yields  a  strain  resolution  of 
about  1 4  ppm. 

The  offset  of  the  horizontal  beams  in  Figure  4.1  is  critical  to  this  d/e  relationship. 

To  examine  the  effect  of  possible  fabrication  flaws  in  this  offset,  we  consider  a  simple 
geometric  construction  in  which  the  horizontal  beams  are  attached  to  the  main  vertical  beam 
by  freely  moving  hinges.  The  resulting  relationship  is: 

(4.2)  d  —  (SL(-)/5)  £  , 

where  d  and  £  are  as  above,  S  is  the  lengdi  of  the  vertical  beam  from  tip  to  midway  between 
the  horizontal  beam  attachment  points  (the  "hinges"),  Lq  is  the  initial  distance  between  the 
anodically  bonded  anchors  (the  distance  labeled  "1mm"  in  Figure  4.1),  and  5  is  the  offset 
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distance  between  hinges.  Taking  S  to  be  535  [im,  Lq  to  be  1  mm,  and  5  as  its  nominal  10 
gm,  we  get: 

(4.3)  d  =  (0.054  gm)  e 

This  geometric  model  therefore  predicts  50%  greater  sensitivity  to  strain  than  does  the  finite 
element  model  It  has  not  taken  flexure  bending  into  account,  but  is  still  useful  as  a  model 
for  examining  the  sensitivity  of  d/ £  to  the  offset  6. 

The  pointer  beam  devices  used  in  the  gamma  irradiations  were  fabricated  at  CSDL 
with  a  0.1  pm  tolerance  for  in-plane  distances.  Examining  equation  4.2,  we  find  that  a  0.1 
pm  change  in  6  results  in  a  1%  change  in  d/£.  A  1%  variation  among  pointer  beam  strain 
sensitivities  would  not  significandy  affect  the  results  we  obtained. 

Three  wafers  were  made:  one  Pyrex  wafer  and  two  Hoya  SD-2  wafers.  The  devices 
were  arrayed  on  glass  wafers,  among  other  MEMS  diagnostic  devices,  as  shown  in  Figure 
4.2.  The  black  squares  indicate  the  locations  of  test  device  arrays.  At  each  of  these 
locations,  one  pointer  beam  is  present,  giving  a  total  of  57  pointer  beams  per  wafer.  The 
glass  used  was  about  780  pm  diick,  and  the  silicon  was  10  pm  thick. 

4.2  Experimental  procedure 
4.2.1  Pre-irradiation  measurements 

To  read  the  strain  gauge,  one  needs  a  microscope.  A  reflected-light  microscope  with 
a  magnification  of  lOOOX  was  used  to  make  measurements.  The  as-fabricated  beam 
positions  for  each  numbered  device  were  measured  and  recorded.  The  average  starting 
deflection  of  the  pointer  beams  is  shown  in  Table  4.1  (positive  deflection  indicates  the 
silicon  is  in  tension). 
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Pyrex 

+0.34 

Hoy  a  B 

+0.14 

HoyaC 

+0.09 

Table  4.1  Pointer  beam  pre-irradiation  positions 

Tooth  shape  and  orientation  showed  considerable  variation  from  one  device  to  the  next.  As 
a  result,  some  degree  of  interpretation  was  involved  in  making  a  deflection  reading.  In  some 
cases,  it  was  clearly  possible  to  make  fractional- tooth  readings.  When  time  allowed,  multiple 
measurements  of  the  same  pointer  beam  were  made,  and  the  results  averaged.  One  pointer 
beam  on  the  Pyrex  wafer  was  found  to  be  at  0.85  teeth.  Phis  pointer  beam  was  discounted 
as  defective,  and  not  used. 
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4.2.2  University  of  Massachusetts  at  Lowell  facility 

A  cobalt-60  gamma  source  at  the  University  of  Massachusetts  at  Lowell  was  used  to 
irradiate  the  pointer  beam  samples.  The  source  consisted  of  a  61  cm  X  61  cm  plaque  of  ^^Co 
strips.  The  strips  were  of  var)ing  activities.  To  irradiate  the  "gamma  cave"  (the  room  in 
which  the  samples  sat),  the  ^'‘to  plaque  would  be  moved  into  a  holder  on  an  aluminum 
window  that  separated  the  source  from  the  samples.  The  gamma  cave  was  a  very  humid 
environment  because  the  concrete  wall  with  the  aluminum  window  also  served,  on  its  other 
side,  as  the  wall  of  a  large  water  tank  in  which  the  ^^Co  source  was  kept.  The  air  inside  was 
sufficiendy  humid  that  water  condensed  on  the  aluminum  window  and  sample  holder. 


Figure  4.3  Pointer  beam  holder 
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4.2.3  Pointer  beam  wafer  holder 


An  aluminum  holder  for  the  samples  was  fasliioned  and  put  into  position  against  the 
aluminum  window  as  shown  in  Figure  4.3.  An  aluminum  sheet  1.6  mm  thick  was  bent  into 
an  L-shape,  then  brads  were  fixed  onto  it  so  tliat  the  wafers  could  be  slid  from  the  top  of  the 
holder  down  to  where  they  were  supported  by  the  brads.  Edges  on  the  brads  prevented  the 
wafers  from  tipping  off  the  holder.  The  MEMS  devices  faced  away  from  the  holder.  The 
wafers,  designated  Pyrex,  Hoya  B,  and  Hoya  C,  were  set  into  the  left,  middle,  and  right 
holder  positions,  respectively. 


Figure  4.4  Placement  of  dosimeters 

4.2.4  Dosimetry 

Prior  to  the  irradiation,  dosimetry  measurements  were  made  by  Lowell  personnel 
using  FWT-70-40M  Optichromic  Dosimeters.'^  ^  Six  dosimeter  packages  were  placed  on  the 
sample  holder,  as  shown  in  Figure  4.4  by  the  dark  rectangles  labeled  as  series  "a".  Each 
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dosimeter  package  held  two  dosimeters,  the  readings  of  which  were  averaged  to  give  one 
reading  per  package.  The  results  are  shown  in  Figure  4.5  with  tlie  package 
centers'  horizontal  positions  measured  with  respect  to  the  center  of  the  middle  wafer.  The 
data  appear  to  follow  a  roughly  gaussian  trend,  and  a  fitted  gaussian  curve  is  shown  in  Figure 
4.5. 


Figure  4.5  Dose  rate  at  the  sample  holder 

A  second  dosimetry  measurement  was  made  at  a  later  date  in  which  four  dosimeters 
were  distributed  vertically  across  the  center  wafer  position  and  two  more  dosimeters  were 
placed  at  the  top  and  bottom  of  the  right  wafer  position.  This  arrangement  is  shown  in 
Figure  4.4  by  tighter  colored  rectangles  labeled  series  "b".  The  measured  dose  rates  were 
found  to  be  constant  to  within  5%  on  the  center  wafer,  and  within  10%  on  the  right  wafer. 
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The  results  of  both  dosimetry  measurements  are  summarized  in  Table  4.2. 


Dosimeter 

Dose  rate,  krads 

Dosimeter 

Dose  rate,  krads 

la 

260 

lb 

720 

2a 

500 

2b 

708 

3a 

672 

3b 

744 

4a 

640 

4b 

744 

5a 

384 

5b 

612 

6a 

164 

6b 

564 

Table  4.2  Dosimeter  readings 

Given  the  placement  of  die  wafers  near  the  vertical  center  of  the  aluminum  window 
(near  the  flat  peak  of  the  assumed  gaussian  vertical  variation),  doses  were  determined  for 
individual  pointer  beams  based  solely  on  their  horizontal  posidons  and  the  fitted  dose  rate 
curve  in  Figure  4.5.  The  average  dose  rate  measured  at  the  right  wafer  position  differed  by 
about  50%  in  the  two  measurement  series,  though  the  center  wafer  position  measurements 
remain  within  7%  of  each  other.  Since  the  wafer  positions  are  not  symmetric  about  the 
gaussian  in  Figure  4.5,  it  is  suspected  tliat  in  the  second  dosimetry  measurement  the  Co-60 
source  plaque  was  placed  in  its  holder  with  an  opposite  orientation  with  respect  to  the  first 
set  of  measurements.  The  ^'‘to  source  was  removed  and  replaced  several  times  over  the 
course  of  the  irradiations  without  verification  of  the  orientation  of  the  source,  so  a  large 
uncertainty  in  the  radiation  dose  received  for  the  left  and  right  wafers  must  be  assumed.  The 
effect  of  this  different  source  placement  is  considered  in  the  discussion  below. 

The  doses  quoted  from  the  dosimeter  measurements  are  water-equivalent  doses  (the 
dose  that  would  have  been  absorbed  by  water  at  the  dosimeter  positions).  To  find  the  doses 
absorbed  by  Pyrex  and  Hoya  SD-2  glasses  for  the  same  gamma  flux,  a  conversion  factor  was 
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calculated.  A  Co-60  gamma  emission  spectrum  of  two  lines  of  equal  intensity,  1.17  MeV  and 
1.33  MeV,  was  assumed.'^^^  Degradation  of  the  emission  spectrum  upon  propagation 
through  the  aluminum  window  and  sample  holder  was  accounted  for  using  mass  attenuation 
coefficients  from  NIST's  website,^  ’  then  the  absorbed  dose  was  calculated  using  mass 
energy-absorption  coefficients,  also  from  NIST."^^  Values  for  absorbed  dose  were 

next  calculated  for  water  and  the  glasses,  then  compared.  Pyrex  and  Hoya  SD-2  absorbed 
doses  were  both  found  to  be  0.89  times  those  measured  by  the  dosimeters  for  water.  It 
turns  out  that  the  thin  aluminum  windows  present  have  a  negligible  effect  for  present 
purposes  on  the  ^^to  emission  spectrum,  and  the  conversion  factor  can  simply  be  taken  as 
the  ratio  of  the  \l^,J p  values  for  the  materials  involved. 

4.2.5  Irradiations  and  measurements 

The  irradiation  of  the  wafers  was  often  interrupted  by  the  removal  of  the  source 
in  order  to  perform  other  simultaneous  experiments  in  the  gamma  cave.  The  wafers  were 
taken  to  CSDL  for  measurement  during  these  interruptions.  The  pointer  beams  on  all  three 
wafers  were  all  measured  after  317,  471,  and  630  hours  of  irradiation.  The  Hoya  C  wafer 
had  its  pointer  beams  additionally  measured  after  184  hours  of  irradiation. 

4.3  Results 

The  results  of  the  pointer  beam  measurements  are  summarized  in  Figure  4.6  and 
Figure  4.7. 
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Figure  4,6  Gamma-induced  strain  in  a  Pvrex  elass  wafer 


A  line  has  been  fitted  to  the  data  using  a  least-squares  technique.  Error  bars  are 
shown  with  ±  1  XlO'*^  strain,  which  is  slightly  better  than  the  14  ppm  resolution  of  the  vernier 
on  the  pointer  beams.  Table  4.3  shows  a  summary  of  the  fits  to  the  data. 


Strain/ rad 

Intercept  (strain) 

Reduced 

Pyrex 

9.02  X  10  ’^ 

2.13  X  10  '’ 

0.65 

Hoya  SD-2 

4.34  X  10  ’" 

-8.51  X  10" 

0.48 

Table  4.3  Pointer  beam  gamma  irradiation  results 

The  different  dosimetry  profile  that  would  result  if  the  ^'^Co  source  were  inserted  into 
its  holder  in  an  orientation  inverted  about  the  vertical  axis  in  Figures  4.3  and  4.5  has  been 
considered.  Dose  values  were  reassigned  to  individual  pointer  beams  based  on  the  curve 
that  results  from  inverting  the  gaussian  curv^e  shown  in  Figure  4.5  about  the  ordinate  axis.  A 
line  was  then  fitted  to  both  the  Pyrex  and  the  Hoya  SD-2  data.  The  resulting  slope  for  the 
Hoya  SD-2  data  was  4.65  X  10  strain  per  rad,  and  1.46  X  10  strain  per  rad  for  the  Pyrex 
data.  This  is  a  7%  difference  for  Hoya  SD-2  and  a  factor  of  6  difference  for  Pyrex  from  the 
values  in  Table  4.3.  These  values  represent  the  extreme  case  that  the  source  sat  for  the 
entire  irradiation  duration  in  the  opposite  orientation  than  previously  thought.  Table  4.4 
summarizes  these  fits,  and  Figures  4.8  and  4.9  show  the  pointer  beam  data  replotted  for  this 
alternate  dosimetry. 


Strain/ rad 

Intercept  (strain) 

Reduced 

Pyrex 

1.46  X  10  ’^ 

3.39  X  10  '’ 

0.66 

Hoya  SD-2 

4.65  X  10  ’" 

-2.35  X  10  '’ 

0.49 

Table  4.4  Pointer  beam  gamma  irradiation  results  (alternate  orientation) 
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Table  4.5  compares  the  fitted  slopes  (strain  per  rad)  for  the  two  extreme  cases  of 


constant  source  orientation  and  presents  the  more  likely  averaged  values. 


Strain/ rad 

Source  orientation  #1 

Strain/ rad 

Source  orientation  #2 

Strain/ rad 

Average 

Pyrex 

9.02  X  10  ” 

1.46  X  10  ” 

5.24  (±  3.78)  X  10  ” 

Hoya  SD-2 

4.34  X  10  ” 

4.65  X  10  ” 

4.50  (±  0.16)  X  10  ” 

Table  4.5  Pointer  beam  gamma  irradiation  averaged  results 

4.4  Discussion 

While  a  trend  of  rising  positive  strain  with  increasing  dose  is  clearly  seen  in  only  the 
Hoya  SD-2  data,  it  is  statistically  shown  to  be  present  in  both  glasses.  The  reduced  chi- 
squared'^'^  values  for  the  fits  indicate  with  high  certainty  that  an  upward  trend  is  present. 

These  results  apparently  contradict  measurements  in  Pyrex  by  Shelby'^  \  At  100 
Mrad,  Shelby  shows  strain  of  approximately  -20  ppm,  whereas  the  present  data  show  a  strain 
in  the  Pyrex  of  about  +3  ppm.  As  Shelby’s  results  have  been  confirmed  by  later 
experiments,'^  '  there  must  be  some  other  element  acting  in  the  pointer  beam  experiment 
to  give  results  so  clearly  different. 

Several  explanations  for  the  Pyrex  results  have  been  proposed.  Firstly,  radiation 
might  be  compacting  the  silicon  pointer  beams  along  with  the  glass  substrate.  In  pure 
silicon,  the  radiolysis  mechanism  involving  Si-O  bonds  is  absent,  and  the  silicon  is  therefore 
not  susceptible  to  the  radiolytic  damage  caused  by  y  rays  in  glass.  Compton  scattered  atomic 
electrons  with  kinetic  energy  equal  to  the  incoming  y  ray  can  be  produced,  however.'^  ^ 
Bombardment  of  silicon  with  1  MeV  electrons  shows  that  expansion  in  silicon  occurs.'^  ^  Not 
only  should  this  expansion  increase  the  differential  strain  with  respect  to  a  compacting  glass 
substrate,  but  its  magnitude  is  negligibly  small  when  to  compared  to  the  expected  Pyrex 
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compaction. 


Secondly,  one  must  consider  diat  the  glass  being  irradiated  has  been  through 
bonding  and  microfabrication  processes  during  the  construction  of  the  pointer  beams.  The 
bonding  process,  called  anodic  bonding,  is  particularly  suspect  as  a  cause  for  changes  in 
material  properties.  During  anodic  bonding,  the  glass  is  heated  to  350-400  °C  while  in 
contact  with  a  silicon  wafer.  Electrodes  on  either  side  of  the  two  wafers  are  set  at  a  potential 
difference  of  500-1000  V.  The  large  electric  field  that  results  causes  sodium  ions  to  drift 
away  from  the  bonding  interface  as  a  layer  of  Si02  forms  at  the  interface  of  die  two 
wafers.'* The  Vickers  hardness  of  Pyrex  has  been  observed  to  increase  two-fold  after 
anodic  bonding.'""  It  is  possible  that  the  chemical  and  microstructural  changes  induced  in 
Pyrex  by  the  anodic  bonding  dramatically  change  the  response  of  the  glass  to  gamma 
irradiation. 

Thirdly,  the  pointer  beams  may  not  be  performing  as  designed.  To  verify  the 
functioning  of  the  pointer  beams,  a  thermal  calibration  apparatus  has  been  designed  and 
built.  An  aluminum  plate  with  a  thermistor,  a  temperature  controller,  and  heaters  attached 
to  it  will  heat  die  glass  wafer  from  underneath  while  a  silicon  circumferential  ring  with  a  glass 
wafer  resting  on  it  act  as  an  oven  for  the  pointer  beams.  Heating  the  pointer  beams  and  the 
wafer  to  which  they  arc  bonded  causes  a  differential  strain  to  develop  between  the  pointer 
beams  and  the  glass  substrate  due  to  differing  coefficients  of  thermal  expansion.  By 
measuring  beam  deflection  and  temperamre,  one  can  calibrate  the  beam  deflection-to-strain 
relationship.  Unfortunately,  this  experiment  is  still  awaiting  a  microscope  powerful  enough 
to  image  the  vernier  on  the  pointer  beams,  but  with  an  objective  lens  able  to  withstand  the 
heat  of  the  experiment. 

Lastly,  effects  such  as  radiation-induced  creep,  interface  debonding,  and  non- 
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dcnsifying  deformation  may  be  possible.  Experiments  involving  heavy  ion  bombardment  of 
silicon,  silica,  and  Pyrex  have  shown  that  non-densifying  radiation-induced  viscous  flow  can 
occur"*  At  high  enough  doses,  the  dimensional  expansion  of  silica  and  Pyrex  (even  in  the 
absence  of  stress)  perpendicular  to  the  irradiating  ions  can  reverse  the  shrinkage  due  to 
saturating  compaction  that  is  evident  at  lower  doses."*  *^’"*  *"*  The  expected  compaction  of  the 
glass  substrate  is  providing  a  stress  that  may  assist  the  radiation-induced  viscous  flow  in  a 
thin  Si02  layer  at  the  anodic  bond,  thus  allowing  the  pointer  beam  anchors  to  slip,  removing 
the  strain  in  the  silicon  pointer  beams.  This  would  not  explain  the  slight  tension  present  in 
the  beams  after  irradiation,  and  the  further  study  must  be  done  to  determine  if  the  gamma 
irradiation  is  of  sufficient  energy  and  dose  to  cause  viscous  flow  under  stress. 
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5  NEUTRON  IRRADIATION  OF  GLASS 


5.1  Overview 

While  the  TRIM  calculations  discussed  in  chapter  3  serve  to  estimate  displacements 
in  fast  and  thermal  neutron  collision  cascades,  they  do  not  indicate  how  many  of  these 
displacements  are  permanent,  nor  how  the  number  of  displacements  affects  the  expansion 
or  compaction  of  a  material  In  order  to  determine  the  effect  of  fast  and  thermal  neutrons 
on  the  density  of  Pyrex  and  Hoya  SD-2  glass,  glass  samples  were  exposed  to  neutron 
fluences  and  the  resulting  change  in  their  density  was  measured  using  a  technique  called 
sink/ float. 

5.2  Experimental  procedure 
5.2.1  Sample  preparation 

Pyrex  and  Hoya  SD-2  wafers  approximately  0.78  mm  thick  were  diced  using  a 
diamond  saw.  Squares  approximately  3  mm  X  3  mm  were  cut  for  use  in  the  sink/ float 
apparatus.  Sink/float  compares  the  density  of  one  material  sample  to  that  of  another. 
Accordingly,  samples  were  needed  both  to  be  irradiated  and  to  be  left  unirradiated  for 
comparison  purposes.  For  each  planned  exposure  level,  samples  were  prepared  in  pairs, 
with  one  to  be  irradiated  and  one  not.  To  minimize  the  probability  of  the  glass  pieces  within 
a  sample  pair  having  an  initial  density  difference,  care  was  taken  to  take  sample  pairs  from 
the  same  wafer,  and  even  to  take  adjacent  pieces  of  glass  as  pairs.  Figure  5.1  shows  an 
example  of  how  the  pairs  were  selected  from  a  Pyrex  wafer. 
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Figure  5.1  Cut  Pyrcx  wafer  and  sample  assignments  (R:  reference) 

In  anticipation  of  the  need  for  distinguishing  between  the  irradiated  samples  and 
their  reference  samples,  the  reference  samples  had  a  corner  cut  off.  While  some 
discok>ration  due  to  irradiation  was  expected,  whether  or  not  this  would  suffice  for 
distinguishing  the  two  samples  during  a  sink/float  measurement  was  not  known.  Figure  5.2 
shows  the  dimensions  of  an  irradiated  sample  and  its  companion  reference  sample. 
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3  mm 

Figure  5.2  Glass  sample  pair 

5.2.2  Neutron  irradiations 

Samples  were  irradiated  by  MIT  nuclear  reactor  personnel  in  a  pneumatic  tube, 
designated  as  IPHl,  inserted  into  the  core  of  the  MIT  reactor.  Estimates  furnished  by  MIT 
personnel  of  the  fast  and  thermal  neutron  fluxes  for  samples  within  IPHl  are  1x10^^ 
n/cm^/s  fast  and  6.7  X  10^^  n/cm^/s  thermal.  Samples  were  held  within  plastic  tubes 
(’’rabbits”)  that  are  2.5  cm  in  diameter  and  about  6.1  cm  long.  The  temperature  during 
irradiation  was  not  monitored,  but  did  not  exceed  the  softening  point  90  ®C)  of  the 
polystyrene  supports  that  held  the  samples  within  the  rabbits,  as  evidenced  by  the  lack  of 
deformation  in  any  of  the  supports  upon  withdrawal  of  the  rabbits  from  IPHl. 

An  irradiation  schedule  was  decided  for  tlie  glass  samples  based  on  work  by 
Paymafl  ^  and  the  TRIM  calculations  discussed  in  chapter  3.  Based  upon  the  experiment 
design  (discussed  below),  it  was  determined  that  a  strain  in  compaction  of  1  X  10  was  the 
greatest  that  would  be  practical  to  measure.  Paymal’s  data  for  Pyrex  compaction  was  for 
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thermal  neutrons  only,  and  predictions  made  based  on  the  low-flucnce  linear  portion  of  his 
strain-to-fluence  cur\^e  show  that  a  6000  second  exposure  would  result  in  approximately  this 
1  X  10"^  strain  value.  Due  to  the  comparatively  low  fast  neutron  flux  and  the  similarity  of 
thermal  and  fast  neutron  displacement  rates  (100  displacements  per  incident  thermal 
neutron,  22  for  a  fast  neutron),  the  effects  of  fast  neutrons  were  expected  to  be  very  small. 

No  radiation  effects  data  were  available  for  Hoya  SD-2  glass,  so  it  was  decided  that  it 
would  be  exposed  to  the  same  flucnccs  as  the  Pyrex.  While  its  lower  boron  content  should 
result  in  fewer  displacements  due  to  thermal  neutron  absorption,  no  predictions  could  be 
made  as  to  how  much  strain  would  result.  The  radiation  schedule  for  the  glass  samples  is 
shown  in  Table  5.1. 


Pyrex 

Hoya 

SD-2 

Irradiation 
Time,  s 

Thermal  Neutron 
Fluence,  n/ cm^ 

Fast  Neutron 
Fluence,  n/ cm^ 

A 

A 

30 

2.0  X  10'" 

3.0  X  10'^ 

B 

B 

60 

4.0  X  10'" 

6.0  X  10'" 

C 

C 

600 

4.0  X  10'^ 

6.0  X  10'’ 

D 

D 

1800 

1.2  X  10"^ 

1.8  X  10'" 

E 

E 

6000 

4.0  X  lO"’ 

6.0  X  10'" 

Table  5.1  Glass  sample  irradiation  schedule 


5.2.3  Sink/float  measurements 

5.2.3.1  Background  on  the  sink/float  technique 

The  ideas  behind  the  sink/ float  technique  for  measurement  of  density  differences  in 
our  glass  samples  are  primarily  drawn  from  Knight^  ^  Shelby^  and  Conners.^  ^  The  method 
has  also  been  used  by  Sato^’  ^  and  Zdaniewsky^  Paymal^  ^  ^  and  again  by  Shelby^  Shelby 
and  Conners  (who  was  a  graduate  student  working  with  Shelby)  report  the  greatest  precision 
for  the  technique,  with  Shelby^  reporting  50  ppm  precision  for  Ap/ p,  Conners  showing 
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~90  ppin  precision,  and  Sato  claiming  300  ppm  precision.  Converting  to  strain  (A// 1),  these 
precisions  would  be  ~  17  ppm,  30  ppm,  and  100  ppm,  respectively. 

To  perform  a  sink/ float  measurement,  a  pair  of  glass  samples  (one  irradiated,  one 
not,  as  described  above)  are  placed  in  a  test  tube  containing  a  heavy  liquid  in  which  they  just 
float.  The  liquid  is  then  heated  by  means  of  a  water  bath.  If  tire  starting  density  of  the  liquid 
is  sufficiently  close  to  that  of  the  glass,  a  few  degrees  of  heating  will  cause  the  liquid  density 
to  grow  smaller  than  that  of  the  glass  (the  thermal  expansion  of  the  glass  is  much  smaller 
than  that  of  the  liquid).  The  glass  samples  will  now  sink,  and  at  different  temperatures  if 
their  densities  are  different.  By  knowing  the  thermal  expansion  coefficient  for  die  Uquid, 
one  can  relate  the  difference  in  sinking  temperatures  to  a  density  difference. 

An  approximate  formula  that  relates  a  difference  in  sinking  temperature  to  a  density 
difference  can  be  derived  as: 


(5.1) 


Ap  _  3Ar  (at; -a;,) 
D  l-3Arat; 


In  Equation  5.1,  AT  is  the  difference  in  sinking  temperature  between  the  irradiated  and 
unirradiated  glass  samples,  a^;  is  the  linear  coefficient  of  thermal  expansion  of  the  glass,  is 
the  linear  coefficient  of  thermal  expansion  of  the  liquid  (a  third  of  the  volumetric  coefficient 
of  thermal  expansion),  and  p  is  the  glass  density.  Since  the  thermal  expansion  of  the  glass  is 
--lO  '^  and  that  of  the  liquid  is  -10  '',  a  good  approximation  can  be  made  simply  by: 

Ap 


(5.2) 


=  3Ar  at 


Assuming  an  uncertainty  of  0.01  °C  in  each  temperature  reading  and  neglecting  the 
uncertainty  in  the  coefficient  of  thermal  expansion  for  the  liquid,  one  obtains  a  predicted 
uncertainty  in  a  strain  calculation  of  about  6  ppm.  This  is  valid  for  small  strains.  As  the 
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strain  rises,  however,  the  uncertainty  in  the  coefficient  of  thermal  expansion  of  the  liquid 
becomes  important,  giving  a  total  of  about  17  ppm  for  a  1  X  10  strain  in  Pyrex  for  an 
assumed  1  %  uncertainty  in  the  liquid's  coefficient  of  thermal  expansion. 


5.2.3. 2  Apparatus  design 

A  schematic  of  the  sink/ float  device  is  shown  in  Figure  5.3.  Shown  are  six  hea\ty 
liquid-filled  38  mm  o.d.  X  200  mm  Pyrex  test  tubes  (round-bottomed,  in  actuality)  suspended 
by  O-rings  from  holes  in  a  plexiglass  lid  that  are  cut  in  a  regular  circular  pattern.  The  tubes 
arc  stoppered  with  Neoprene®  rubber,  a  material  resistant  to  degradation  through  contact 
with  vapor  from  the  organic  liquids  in  the  test  tube.  The  lid  rests  on  a  cylindrical  25  cm  o.d. 
X  25  cm  Pyrex  jar.  A  2.5  cm  o.d.  X  25  cm  (only  bottom  12.5  cm  heated)  500  watt  immersion 
heater  is  suspended  from  a  hole  in  the  center  of  the  lid.  The  heater  is  attached  to  a  Variac 
transformer.  The  jar  is  filled  with  water  to  above  the  fluid  level  in  the  test  tubes,  and  a  7.6 
cm  stir  bar  is  at  the  bottom  of  the  jar,  for  use  with  a  heavy-duty  stir  plate.  A  horizontal  line 
is  drawn  about  5  cm  up  from  the  test  tube  bottoms. 

A  platinum  resistance  temperature  detector  probe  was  inserted  through  a  rubber 
stopper  and  immersed  as  shown  in  the  liquid  of  one  of  the  test  tubes.  The  probe  was 
attached  to  a  precision  thermometer  with  0.01  readout  and  0.05  accuracy  at  the 
temperatures  of  interest. 
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Figure  5.3  Sink/Float  apparatus  schematic 


5.2.3.3  Set-up  of  sink/ float 

The  glass  samples  were  still  radioactive  when  the  experiment  began.  Only  Pyrex 
sample  E  (see  Table  5.1)  was  above  background  level,  but  all  of  the  Hoya  SD-2  samples 
were  measurably  active,  most  notably  with  ’^^Sb.  It  is  not  certain  whether  the  antimony  is 
intrinsic  to  the  glass  or  if  it  somehow  came  to  be  present  during  handling. 

All  samples  came  back  tinted  a  brown  color,  evidence  of  trapped  electronic  defects. 
The  color  of  the  Pyrex  samples  varied  continuously  and  conspicuously  from  light  to  medium 
brown  as  the  sample  exposure  rose.  In  the  case  of  the  Hoya  SD-2,  the  coloration  seemed  to 
be  cjuite  ecjual  among  the  three  samples  (C,  D,  E)  receiving  the  highest  doses,  which  points 
to  a  possible  saturation  of  the  electronic  defects.  These  observations  are  just  cursory 
impressions,  however,  and  no  quantification  of  this  effect  was  made  using  optical  absorption 
spectrometry  or  other  techniques. 

Samples  were  cleaned  by  rinsing  them  with  acetone.  During  the  cleaning  of  tlie 
Hoya  SD-2  samples,  problems  with  air  currents  in  the  vent  hood  led  to  a  complete  mix-up 
of  samples.  The  relative  activities,  however,  helped  to  sort  them  out.  Their  activities  had 
been  individually  measured  upon  their  arrival  at  CSDL,  and  monotonically  rose  with 
exposure  time. 

The  heavy  hquids  used  were  tetrabromoethane  (also  known  as  acetylene 
tetrabromide)  and  bromonaphthalene.  The  first  was  suspected  of  being  quite  hazardous,  so 
all  work  was  performed  in  a  vent  hood  while  wearing  goggles,  a  respirator  with  organic 

vapor  filter  cartridges,  and  Viton  gloves. 

The  tetrabromoethane  has  a  density  of  2.96  g/ cm  ,  and  the  bromonaphthalene  has  a 
density  of  1.48  g/evn.  Mixtures  of  these  liquids  were  made  in  a  IL  beaker.  Initial  liquid 
proportions  were  based  on  the  manufacturers  densities  cited  here  for  tlie  liquids  and  the 
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glasses  (pp,,,,  2.23  g/cm'\  p|,,,sD-2  =  2.60  g/cm').  All  of  the  sample  pairs  were  individually 

dropped  in  the  mixture  to  verify  that  they  all  floated,  and  the  liquid  was  adjusted  by  adding 
milliliter  quantities  of  die  lighter  or  denser  liquid.  The  densest  sample  was  identified  and 
then  just  made  to  float. 

5.2.3.4  Sink/ float  execution 
5.2.3.4.1  Pyrex 

Measurement  runs  were  first  performed  on  Pyrex.  The  apparams  had  equilibrated 
overnight  at  about  22  °C.  After  getting  a  feel  for  the  temperature  response  to  heating,  the 
transformer  output  was  set  to  54  volts,  with  a  resulting  observed  current  of  0.4  amperes,  or 
22  watts.  Tliis  gave  a  0.11  ±  0.02  °C  per  minute  liquid  heating  rate.  The  liquid  temperaUire 
rose  steadily  and  predictably  —  a  given  0.01  °C  increment  is  traversed  in  about  6  seconds. 

The  system  showed  a  response  time  of  about  10  minutes  to  a  change  in  heater  power.  A 
high  stir  setting  was  used,  with  the  presence  of  the  mbes  preventing  a  vortex  from  forming 
(with  the  tubes  removed  at  this  setting,  the  vortex  extends  to  the  very  bottom  of  the  jar). 

A  sample  was  declared  as  "sunk"  when  it  crossed  the  line  drawn  on  the  test  tubes. 
Samples  tended  to  stick  to  the  surface  of  the  hquid  (and,  less  often,  to  each  other),  and  only 
very  vigorous  shaking  freed  them  up.  Samples  took  several  minutes  to  sink  once  they  had 
left  the  surface,  sometimes  hovering  at  a  fixed  depth  before  resuming  the  fall. 

Sinking  of  samples  usually  occurred  in  the  center  of  the  tube,  and  convective 
currents  within  tire  test  tubes  appeared  to  be  present.  After  sinking  and  hitting  the  bottom 
of  the  test  tube,  Pyrex  sample  E  (the  densest,  most  irradiated)  would  move  to  the  side  of  the 
tube  and  rise  again  about  2/3  of  the  way  to  the  top.  It  would  repeat  this  cycle  three  times, 
with  each  subsequent  rise  height  being  smaller  than  the  one  before,  then  stick  at  the  bottom. 
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In  all  cases,  the  recorded  sinking  temperature  recorded  was  that  of  a  sample's  first  crossing 
of  the  reference  line. 

Viewing  of  tlic  samples  was  clear  and  easy.  This  was  greatly  helped  by  the 
magnifying  effect  of  die  water-filled  cylindrical  jar  and  the  heavy  liquid-filled  test  tubes.  The 
irradiation-induced  darkening  of  samples  C,  D,  and  E  would  have  been  sufficient  to  easily 
distinguish  between  irradiated  and  unirradiated  glass  samples,  but  for  samples  A  and  B,  the 
cut  corner  on  the  unirradiated  sample  was  necessary.  A  typical  falling  position  was  edge-on, 
so  the  cut  corner  probably  made  minimal  difference  in  the  drag  characteristics  of  the  sample. 

"First  fall"  (of  sample  E)  occurred  between  28.15  °C  and  28.65  °C.  It  took  about 
one  hour  to  heat  to  this  point.  Sinking  temperatures  for  the  unirradiated  glass  occurred  in 
the  31.89-32.39  ""C  range.  The  temperatures  are  those  read  in  one  test  tube  -  the  others 
probably  do  not  share  this  temperature,  but  it  is  the  difference  in  sinking  temperatures  that  it 
is  desired  to  know.  Sinking  temperatures  varied  for  a  given  sample  in  a  given  tube  to  within 

about  a  0.1  °C  standard  deviation. 

On  subsequent  runs,  the  water  bath  was  emptied  by  siphoning  then  refilled  with  tap 
water  (colder  than  room  temperature).  Rapid  heating  was  applied  to  about  24  °C,  then  the 
22  watts  setting  was  resumed.  Given  the  observed  --  10  minutes  needed  to  achieve  steady- 
state  temperature  rise,  this  left  plenty  of  time  before  "first  fall". 

Runs  were  made  in  which  the  starting  temperature  was  high  and  all  samples  were 
firmly  at  the  bottom  of  tlic  test  tubes  at  the  start.  By  stirring  without  heating,  the  liquid 
temperature  cooled.  The  cooling  rate  obtained  this  way  was  rmpractically  low,  so  cool  water 
was  added.  This  process  was  not  well  controlled,  and  resulted  in  results  that  do  not  at  all 
agree  with  the  heating  runs. 
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Runs  were  made  in  which  tlie  tube  shaking  was  varied  as  to  degree  —  there  was  some 
concern  that  the  jiggling  of  test  tubes  was  a  contributor  to  currents  within  the  heavy  liquids. 

A  run  made  with  no  shaking  clearly  delayed  the  sinking  of  samples  due  to  the  liquid  surface- 
tension  attraction.  Five  of  the  Pyrex  runs  were  done  with  an  "only  as  necessary”  test  tube 
shaking  policy;  tubes  were  jiggled  about  every  minute  until  the  samples  started  to  sink. 
Once  both  samples  in  a  tube  were  detached  from  the  liquid  surface,  no  further  shaking  took 
place.  A  run  was  done  with  constant  shaking  of  die  test  tubes  throughout  the  run,  and  the 
results  are  in  agreement  with  five  other  runs. 

5.2.3.4.1  Hoya  SD-2 

The  Hoya  SD-2  measurements  proceeded  as  above,  with  two  differences  to  note. 
First,  the  attraction  of  the  surface  for  the  samples  seemed  somewhat  less  than  for  the  Pyrex 
case.  Second,  the  liquid  density  required  significantly  more  adjustment  than  that  of  the 
Pyrex  -  based  on  the  published  liquid  densities  and  the  final  proportions  in  the  mixture,  die 
Hoya  SD-2  glass  appears  to  have  a  density  of  2.62  g/ cm^.  About  62  ml  extra 
tetrabromo ethane  was  needed  to  get  all  the  samples  to  float,  then  5  ml  extra 
bromonaphthalene  was  added  ("extra”  referring  to  a  planned  density  of  2.60  g/ cm'^  and  780 
ml  total  fluid).  As  a  result  of  more  careful  mixing,  the  liquid  mixture  was  much  closer  to  the 
"first  fall"  point  than  for  the  Pyrex.  In  the  first  run,  starting  from  22.14  °C,  first  fall 
occurred  at  22.53  °C.  There  was  concern  that  this  sinking  temperature  so  close  to  the 
starting  point  wasn't  allowing  steady  state  to  be  reached,  so  chilled  water  was  added  to  the 
water  bath  on  subsequent  runs.  Start  temperatures  for  the  heavy  liquid  were  then  close  to  20 
°C.  The  results  from  the  first  run  resemble  those  of  subsequent  runs,  so  have  been  kept. 
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All  runs  were  performed  with  the  above-described  "only  as  necessary"  tube  shaking 


approach. 

5.3  Results  and  discussion 

The  results  for  six  measurements  each  of  Pyrex  and  Hoya  SD-2  are  shown  in  Figure 

5.4  and  Figure  5.5,  respectively.  The  error  bars  shown  represent  the  standard  deviation  of 
the  six  measurements  divided  by  the  square  root  of  six,  then  averaged  across  all  five  samples. 
This  is  the  uncertainty  for  the  mean  value  of  the  six  measurements.^  It  is  important  to 
distinguish  between  this  scatter,  which  is  indicative  of  the  repeatability  and  precision  of  the 
sink/float  method,  and  the  acci/rag  of  the  conversion  of  a  measured  AT  to  a  strain.  An 
inaccuracy  in  the  thermal  expansivity  of  the  heavy  liquids  used  will  add  to  tlie  inaccuracy  of 
the  strain  values,  but  will  not  contribute  to  scatter.  Scatter  in  the  points  is  due  solely  to 
uncertainty  in  temperature  readings,  convective  currents,  and  delays  in  sinking  due  to  liquid 
surface  attraction.  In  the  ideal  situation  (no  current-induced  inaccuracies,  no  surface- 
sticking),  and  with  our  0.01  ""C  temperature  control,  we  expect  6  ppm  scatter  in  the  data. 
Table  5.2  summarises  the  observed  scatter. 


Pyrex 

Hoya  SD-2 

Average  a 

55  ppm 

16  ppm 

Average  _5_ 

V6 

22  ppm 

6  ppm 

Table  5.2  Average  scatter  a  of  six  measurements  and  resultant  uncertainty  of  the  mean 
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Figure  5.5  Neutron  irradiation-induced  strain  in  Hoya  SD-2  as  measured  by  sink/ float 


The  Pyrex  scatter  has  an  order  of  magnitude  greater  scatter  than  the  ideal  case,  yet  is 
similar  to  the  results  attained  by  Shelby  and  Conners.  The  Hoya  SD-2  data  is  more  tighdy 
grouped.  This  may  be  related  to  the  above  observ^ation  that  these  samples  seemed  to  stick 
less  to  the  liquid's  surface.  The  average  of  six  points  attains  the  precision  of  the  ideal 
sensitivity  of  a  single  measurement. 

Additional  uncertainty  (but  not  scatter)  is  introduced  by  the  uncertainty  of  the 
thermal  expansivity  of  the  liquid.  The  data  analysis  here  has  been  done  using  a  value 
measured  by  Conners^  ’.  Conners  used  sink/ float  to  detect  gamma-induced  density  changes 
in  glasses  with  densities  from  2.13  g/ cm^  to  2.52  g/ cm^.  Her  liquid  mixtures  were  all 
composed  of  the  same  two  chemicals  used  here,  tetrabromoethane  and  bromonaphthalene. 
She  observed  that  the  volumetric  thermal  expansion  for  all  of  her  liquids  averaged  to  840 
ppm/°C,  with  a  standard  deviation  of  8  ppm/°C,  and  this  was  the  value  she  then  applied  to 
all  of  her  data  analysis.  A  1%  uncertainty  in  the  value  for  the  thermal  expansivity  of  the 
liquid  introduces  a  1  %  uncertainty  into  the  level  of  strain  being  measured. 

The  linear  fit  shown  gives  a  reduced  for  four  degrees  of  freedom  of  1.01  for  the 
Pyrex  and  1.58  for  the  Hoya  SD-2,  indicating  reasonable  fits  to  the  data.^  ”  The  fit  was  made 
to  tlie  averaged  points  using  the  values  in  the  second  row  of  Table  5.2  as  uncertainties,  and 
included  the  zero  point,  to  which  an  uncertainty  equal  with  that  of  the  other  fitted  points 
was  assigned.  Figures  5.6  and  5.7  show  an  enlarged  view  of  the  low-irradiation  time  data, 
and  Figure  5.8  shows  a  comparison  between  the  Pyrex  and  the  Hoya  SD-2  results.  Table  5.3 
presents  a  summary  of  the  linear  fits. 


79 


4.0E-05 


^  n  rO  04  fO 

^ 


o 

o 

3 


»n  lT)  o 


D  :3 
VO  VO 


□  O  <  O  X  4- 


<  i£! 


ureJJS  z-as 


0  100  200  300  400  500  600  700 

Irradiation  Time,  seconds 

Figure  5.7  Neutron  irradiation-induced  strain  in  Hoya  SD-2  as  measured  by  sink/ float  technique 
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Figure  5,8  Neutron  irradiation-induced  strain  in  Pyrex  and  Hoya  SD-2  as  measured  by  sink/ float 


Slope,  strain/ s 

Intercept 

Reduced 

Pyrex 

-1.85x10’ 

1.19x10'^ 

1.01 

Hoya  SD-2 

-6.82x10^ 

5.68x10" 

1.58 

Table  5.3  Linear  fitting  results  for  sink/ float  data 


Using  the  slopes  from  the  fits  shown  in  Table  5.3  and  combining  this  with  the  fast 
and  thermal  neutron  collision  cascade  information  from  Chapter  3  and  the  glass 
compositions  and  densities  presented  in  Chapter  2,  we  obtain  the  values  shown  in  Table  5.4 
for  thermal  neutron-  and  fast  neutron-induced  strain.  Neutron  incidence  into  a  thickness  of 
0.78  mm  is  assumed. 


strain/ dpa 

strain/ n/ cm^ 
thermal 

Strain/ n/ cm^ 
fast  (1  MeV) 

Pyrex 

-1.49 

-2.77x1 0'^" 

-6.07x10"' 

Hoya  SD-2 

-0.19 

-1.01x10'^' 

-7.89x10'"" 

Table  5.4  Strain  in  glass  related  to  neutron  fluence  and  displacements  per  atom 
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6  DISCUSSION 


The  results  for  thermal  neutron  damage  Pyrex  agree  closely  with  those  obtained  by 
Paymaf’  tliough  there  is  some  question  as  to  whether  the  results  are  direcdy  comparable, 
as  Paymal  does  not  explicitly  state  what  his  sample  dimensions  are.  He  does  give  several  of 
clues,  but  they  seem  to  be  contradictory.  He  mentions**^  a  Pyrex  sample  size  of  5  mm  X  5 
mm  X  5  mm  in  connection  with  post-irradiation  thermal  annealing  studies,  also  using  the 
sink/float  method,  but  does  not  connect  this  with  the  fluence  vs.  density  change 
measurements.  As  to  the  latter,  he  claims''  ’  without  justification  that  one  can  estimate  the 
number  of  thermal  neutrons  absorbed  by  the  Pyrex  according  to  the  relationship:  n=3D, 
where  n  is  the  number  of  boron-absorbed  neutrons  per  cm\  and  D  is  the  number  of 
incident  thermal  neutrons  per  cm^.  My  calculation  based  on  the  boron  density  in  Pyrex  and 
its  cross-section  for  thermal  neutron  capture  yields  a  relationship  of  n=2.43D  in  the  limit  of 
a  very  thin  sample.  For  a  5  mm  thick  sample,  I  calculate  that  the  relationship  should  be 
n=l  .41D.  In  the  case  of  our  3  mm  x  3mm  x  0.78  mm  samples,  an  average  over  the  three 
dimensions  gives  a  relationship  of  n=1.89D.  We  might  therefore  expect  to  have  1.3  more 
displacements  per  volume  than  Paymal,  and  about  that  much  more  strain.  AU  of  the 
analyses  in  Chapter  5  assume  neutrons  incident  in  the  thin  (0.78  mm)  direction,  and 
therefore  a  uniform  displacement  concentration  for  the  glass.  For  a  thick  material 
unidirectionally  irradiated,  the  ~4  mm  mean  free  path  of  thermal  neutrons  in  Pyrex  (~16 
mm  in  Hoya  SD-2)  can  lead  to  a  strain  gradient  in  the  material. 

Figure  6.1  shows  a  comparison  of  the  sink/float  data  from  the  present  study  to 
Paymal's  data.  As  fast  neutrons  were  present  in  the  current  data,  their  effect  was  accounted 
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Thermal  Neutron  Fluence,  n/cm 
gure  6.1  Comparison  of  thermal  neutron-induced  strain  results 


for  and  removed  by  means  of  ratioing  the  total  strain  by  the  TRIM  numbers  for  fast  and 
thermal  neutron  displacements  from  Chapter  3.  The  calculated  contribution  of  fast 
neutrons  is  only  0.3%  of  the  total  strain.  The  differences  between  the  current  results  and 
Paymal's  are  no  greater  than  those  between  the  results  of  Shelby  and  Conners  discussed  in 
Chapter  2,  and  may  be  related  to  compositional  differences  or  thickness  effects. 

With  such  good  agreement  for  our  neutron  damage  results  in  Pyrex  with  those  of 
previous  studies,  it  can  be  confidently  said  that  Hoya  SD-2  glass  has  greater  resistance  to 
compaction  due  to  atomic  displacement  damage  as  caused  by  fast  and  thermal  neutrons  than 
Pyrex.  The  numbers  shown  in  Table  5.4  clearly  demonstrate  this,  with  Hoya  SD-2 
compacting  almost  an  order  of  magnitude  less  for  a  given  1  MeV  neutron  fluence. 

It  is  clear  that  the  gamma  radiation  effect  on  pre-anodic  bonding  Hoya  SD-2  has  not 
been  established  by  our  results  in  Chapter  4.  At  the  same  time,  we  may  have  evidence  of 
interfacial,  creep,  and/or  viscous  flow  phenomena  in  the  MEMS  devices  used  in  the  gamma 
experiments.  We  may  also  have  evidence  for  the  anodic  bonding  process  fundamentally 
changing  the  radiation  response  of  Pyrex  and  Hoya  SD-2. 

A  comparison  of  the  Pyrex  pointer  beam  results  with  previous  sink/ float  gamma 
irradiadon  studies  (see  Figure  2.4)  is  shown  in  Figure  6.2.  One  might  conclude  that  the 
pointer  beam  results,  wliile  not  showing  the  expected  clear  trend  of  compaction,  do  not 
really  disagree  with  these  other  studies.  Sato’s  data  show  expansion  at  these  doses,  and 
Shelby's  error  bars  almost  overlap  tliose  of  the  pointer  beams.  Shelby  s  data,  however, 
follow  a  very  clean  cur\-e  towards  compaction  that  is  confirmed  by  his  measurements  at 
higher  doses.  At  higher  doses,  Sato's  data  merge  with  Shelby's,  and  with  the  relative 
uncertainty  shown  only  one  of  his  measurements  differs  much  from  Shelby’s.  In  addition,  a 


87 


statistical  analysis  shows  with  high  degree  of  certainty  of  the  presence  of  an  upward  trend  in 
the  pointer  beam  data,  as  discussed  in  Chapter  3. 

A  valid  objection  may  be  raised  as  to  die  validity  of  the  TRIM  calculations  discussed 
in  Chapter  3  for  our  glasses  on  the  grounds  that  atomic  bonding  information  was  left  out  of 
the  glass  model.  Only  relative  magnitudes  of  displacements  were  used,  however,  as  the  data 
were  used  only  to  convert  thermal  neutron  damage  to  fast  neutron  damage  within  a  single 
glass  type.  There  is  some  support  to  be  found  for  the  validity  of  such  comparisons  to  be 
found  in  the  following  observations:  1)  the  fast  neutrons  certainly  did  not  do  significantly 
n7ore  damage  in  comparison  to  thermal  neutrons  than  predicted  by  TRIM,  given  the  close 
agreement  of  our  neutron  compaction  results  in  Pyrex  and  those  of  Paymal  and  2)  the 
conversion  made  in  section  2.3  of  Pyrex  thermal  neutron  strain  to  fast  neutron  strain  using 
TRIM  results  yielded  strain  values  of  Pyrex  vs.  vitreous  silica  that  might  have  been  guessed 
from  the  gamma  irradiation  results  shown  in  Figure  2.6.  It  is  important  to  also  note  that  no 
conclusions  have  been  drawn  directly  from  TRIM  calculations;  no  predictions  have  been 
made  of  strain  based  solely  on  a  dpa  number,  but  have  been  made  for  fast  neutrons  based 
on  thermal  neutron  data  that  have  been  multiplied  by  fast/ thermal  neutron  dpa  numbers 
from  TRIM. 

It  is  important  to  realize  that  the  radiation  doses  absorbed  by  samples  in  the  present 
study  are  well  below  the  saturation  values  shown  for  thermal  neutrons  in  Figure  2.7,  for  fast 
neutrons  in  Figure  2.8,  and  for  gamma  irradiation  in  Figure  2.4  (saturation  of  strain  is  not 
shown  here,  but  a  roU-over  of  the  data  starts  to  happen  at  about  5  X  lO’  rad).  The 
fundamental  relationship  of  strain-to-dose  is  certainly  not  linear.  For  die  relatively  small 
exposures  in  our  experiments,  however,  a  line  seems  as  appropriate  as  anything  else  for 
extrapolating  to  intermediate  or  smaller  doses  of  radiation. 
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7.  CONCLUSION 


7.1  Summary  of  principal  findings 

Density  measurements  for  neutron-irradiated  Pyrex  and  Hoya  SD-2  glass  samples 
were  carried  out.  In  the  analysis  of  the  data,  the  fast  neutron-induced  compaction  was 
separated  from  the  primarily  thermal  neutron-induced  damage  by  ratioing  the  measured 
compaction  by  the  average  number  of  displacements  in  a  fast  neutron  collision  cascade 
relative  to  the  number  of  displacements  for  a  thermal  neutron  absorption  event.  Monte- 
Carlo  calculations  of  the  number  of  displacements  in  each  ty'pe  of  collision  cascade  were 
carried  out  using  TRIM.  Table  7.1  summarizes  the  findings  for  the  relationships  between 
strain  and  neutron  flucncc,  as  well  as  those  between  strain  and  displacements  per  atom  (dpa) 
for  both  glasses.  Figures  7.1  and  7.2  represent  in  graphical  form  the  values  in  Table  7.1. 


Strain/ (n/ cm^) 
thermal 

Strain/ (n/cm^) 
fast  (1  MeV) 

Strain/ dpa 

Pyrex 

-2.77x10'^'’ 

-6.07x10"' 

-1.49 

Hoya  SD-2 

-1. 01X10'^' 

-7.89x10"" 

-0.19 

Table  7.1  Neutron  irradiation-induced  strains 

Silicon  pointer  beam  strain  gauges  bonded  to  glass  wafers  were  irradiated  with 
gamma  rays.  The  differential  strain  between  the  irradiated  Pyrex  wafer  and  pointer  beams 
was  found  to  be  much  smaller  than  was  predicted  based  upon  previous  gamma  irradiation 
studies  of  Pyrex.  It  is  thought  that  tliis  surprising  result  has  its  origin  in  radiation-induced 
creep  at  the  silicon-glass  bond  layer,  and/ or  in  changes  in  the  glass  radiation  response 
induced  by  the  anodic  bonding  process  used  in  the  manufacture  of  the  pointer  beams. 
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Figure  1  SI  Fast  neutron-induced  strain 


The  relationships  between  induced  differential  stein  and  gamma  dose  for  both  glasses  are 
listed  in  Table  7.2. 


Strain/ rad 

Pyrex 

5.24  (±  3.78)  X  10-'^ 

Hoya  SD-2 

4.50  (±0.16)  X  lO  ’" 

Table  7.2  Pointer  beam  gamma  irradiation  results 
7.2  Importance  of  findings  to  MEMS 

At  the  beginning  of  this  study,  our  goal  was  to  identify  the  radiation-induced 
dimensional  changes  in  Pyrex  and  Hoya  SD-2.  Given  a  relationship  of  stein  to  dose  (such 
as  those  listed  in  Tables  1.2  and  1.3),  the  hope  was  that  one  could  determine  the  induced 
strain  in  a  glass  substrate,  as  well  as  the  resulting  stress  in  the  silicon  devices  bonded  to  it. 
With  a  well  defined  radiation  environment  (and  a  knowledge  of  the  silicon  response  to 
radiation),  the  suitability  of  Pyrex  and  Hoya  SD-2  as  substrates  for  a  particular  application 
could  be  determined. 

Our  experiments  and  calculations  produced  values  relating  strain  to  thermal  and  fast 
neutron  fiuence  for  both  glasses.  The  relationship  between  stein  and  gamma  dose  has  been 
investigated  previously  for  Pyrex.  Our  gamma  irradiation  experiment,  aimed  at  producing  a 
similar  relationship  for  Hoya  SD-2,  indicates,  however,  that  the  radiation  effect  on  anodically 
bonded  MEMS  devices  is  not  simply  the  sum  of  the  parts  (tlie  ’’parts”  being  the  separate 
radiation-induced  dimensional  changes  on  the  glass  and  on  the  silicon).  Interfacial  effects  or 
changes  made  in  the  glass  chemical  structure  by  the  anodic  bonding  process  evidently  play  a 
role  in  the  overall  MEMS  radiation  response.  The  relationships  between  dose  and 
dimensional  change  established  here  are  not  only  useful  for  predicting  strains  in  the  separate 
materials,  but  provide  data  necessan^  to  a  better  understanding  of  anodicatly-bonded  MEMS 
devices. 


93 


7.3  Suggestions  for  future  work 

The  pointer  beam  gamma  results  highlight  the  importance  of  studying  the  anodically 
bonded  silicon-glass  system  in  order  to  define  the  roles  of  individual  system  pieces  in  the 
overall  radiation  response.  Irradiation-induced  creep  and  viscous  flow  of  the  bulk  materials 
and  of  the  anodic  bond  interface  layer  should  be  studied.  One  possible  method  might  be  a 
variant  of  the  wafer  cur\"ature  method  used  in  an  ion  irradiation  study  of  vitreous  silica. 
Another  approach  is  fast  and  thermal  neutron  irradiation  of  pointer  beam  devices.  By 
comparing  the  sink/ float  results  with  pointer  beam  results,  one  could  compare  the  effects  of 
radiation  in  as-supplied  glass  to  its  effects  in  the  bonded  MEMS  device  and  substrate. 
Similarly,  glass  samples  should  be  exposed  to  gamma  radiation  and  then  measured  with  the 
sink/ float  technique. 

Sink/ float  experiments  should  be  performed  on  irradiated  (with  gamma  rays  and 
neutrons)  glass  that  has  been  through  the  anodic  bonding  process.  It  is  possible  to 
anodically  bond  a  glass  wafer  to  a  silicon  wafer  and  then  free  them  from  each  other.  In  this 
fashion,  it  would  be  possible  to  measure  the  irradiadon  effects  on  bulk  glass  after  bonding, 
and  isolate  these  effects  from  those  at  the  bond  interface. 

In  addition  to  those  already  carried  out  on  thermal  neutron  irradiated  glasses, 
experiments  involving  fast  neutron  exposures  should  be  carried  out  on  both  glasses  using 
sink/float  measurements  and  pointer  beams.  The  predictions  presented  here  for  fast 
neutron-induced  glass  compaction  rely  on  TRIM  calculations  coupled  with  experimental 
data  for  thermal  neutron  irradiations.  The  results  thus  obtained  should  be  confirmed  by 
direct  experiment. 

Itinally,  the  thermal  calibration  of  the  pointer  beams  mentioned  in  Chapter  4  should 
be  performed  to  confirm  the  proper  function  of  the  pointer  beams. 
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APPENDIX  A 


Pointer  Beam  Strain  Measurements 


Pyrex 


Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

6.30E+07 

-1.79E-06 

1.18E+08 

7.14E-06 

1.53E+08 

-5.36E-06 

2.11  E+08 

-1.43E-05 

6.30E+07 

1.43E-05 

1.18E+08 

-7.14E-06 

1.53E+08 

4.07E-05 

2.11E+08 

3.57E-06 

6.30E+07 

7.14E-06 

1.18E+08 

7.14E-06 

1.58E+08 

-3.57E-06 

2.11E+08 

-7.14E-06 

7.29E+07 

5.36E-06 

1.18E+08 

1.43E-05 

1.58E+08 

-3.57E-06 

2.11  E+08 

-3.57E-06 

7.29E+07 

1.43E-05 

1.18E-I-08 

1.25E-05 

1.58E+08 

O.OOE+OO 

2.12E+08 

1.07E-05 

7.29E+07 

1.43E-05 

1.18E+08 

-1.07E-05 

1.58E+08 

O.OOE+OO 

2.12E+08 

1.43E-05 

7.29E+07 

7.14E-06 

1.18E+08 

-7.14E-06 

1.58E+08 

O.OOE+OO 

2.12E+08 

7.14E-06 

7.29E-f07 

1.43E-05 

1.24E+08 

-3.57E-06 

1.58E+08 

1.25E-05 

2.12E+08 

1.07E-05 

8J6E+07 

-3.57E-06 

1.24E+08 

1.79E-06 

l,58E+08 

7.14E-06 

2.12E+08 

7.14E-06 

8.36E+07 

1.79E-06 

1.24E+08 

-3.57E-06 

1.58E+08 

3.57E-06 

2.12E+08 

1.61E-05 

8.36E-f07 

-3.57E^06 

1.24E+08 

1.07E-05 

1.58E+08 

-3.57E-06 

2.12E+08 

7.14E-06 

8.36E+07 

1.07E-05 

1.24E+08 

^7.14E  06 

1.66E+08 

3.57E-06 

2.12E+08 

1.43E-05 

8.36E+07 

O.OOE+00 

E24E+08 

1.07E-05 

1.66E+08 

1.61E-05 

2.12E+08 

3.57E-06 

8.36E+07 

3.57E-06 

1.24E+08 

O.OOE+00 

1.66E+08 

7.14E-06 

2.27E+08 

-3.57E-06 

8.36E'f07 

O.OOE+00 

1.25E+08 

1.25E-05 

1.66E+08 

1.79E-05 

2.27E+08 

^5.36E-06 

9.36E+07 

5.36E^06 

1.25E+08 

1.43E-05 

1.66E+08 

O.OOE+OO 

2.27E+08 

-7.14E-06 

9.36E-f07 

3.57E-06 

E25E+08 

1.43E-05 

1.66E+08 

1.43E-05 

2.35E+08 

1.07E-05 

9.36E+07 

1.43E-05 

1.30E+08 

-7.14E-06 

1.66E+08 

7.14E-06 

2.35E+08 

1.79E-06 

9.48E+07 

-8.93E-06 

1.30E+08 

1.07E-05 

1.76E+08 

-3.57E^06 

2.35E+08 

7.14E-06 

9.48E+07 

-3.57E-06 

1.30E+08 

7.14E-06 

1.76E+08 

-5.36E-06 

2.35E+08 

-3.57E-06 

9.48E+07 

-3.57E-06 

1.30E+08 

-7.14E-06 

1.76E+08 

O.OOE+OO 

2.35E+08 

1.43E-05 

9.48E+07 

1.07E-05 

1.30E+08 

O.OOE+OO 

1.76E+08 

-7.14E-06 

2.35E+08 

1.79E-05 

9.48E-f07 

3.57E-06 

1.30E+08 

-7.14E-06 

1.76E+08 

7.14E-06  ^ 

2.35E+08 

1.96E-05 

9.48E+07 

O.OOE+00 

1.30E+08 

-3.57E-06 

1.76E+08 

7.14E-06 

2.35E+08 

-7.14E-06 

9.48E+07 

7.14E-06 

i.41E+08 

-7,14E-06 

1.76E+08 

1.25E-05 

2.35E+08 

1.43E-05 

9.48E+07 

3.57E-06 

1.41E+08 

O.OOE+OO 

l,76E+08 

-7.14E-06 

2.59E+08 

1.07E^05 

9.48E+07 

O.OOE+00 

1.41  E+08 

3.57E-06 

1.76E+08 

-3.57E-06 

2.59E+08 

1.43E-05 

1.06E+08 

3.57E~06 

1.41E+08 

3.57E-06 

1.88E+08 

7.14E-06 

2.59E+08 

1.43E-05 

1.06E+08 

^7.14E-06 

1.41E+08 

O.OOE+OO 

1.88E+08 

3.57E-06 

2.59E+08 

7.14E^06 

1.06E+08 

7.14E-06 

1.41E+08 

3.57E-06 

l,88E+08 

1.07E-05 

2.59E+08 

3.57E^06 

1.06E+08 

7.14E-06 

1.41E+08 

3.57E-06 

1.88E+08 

1.43E-05 

2.59E+08 

1.43E-05 

1.06E+08 

7.14E-06 

1.41E+08 

3.57E-06 

1.88E+08 

-7.14E-06 

2.59E+08 

7.14E-06 

1.06E+08 

1.25E-05 

1.42E+08 

-1.79E-06 

1.88E+08 

1.43E-05 

2.82E+08 

1.25E~05 

E06E+08 

7.14E-06 

1.42E+08 

-1.43E-05 

1.88E+08 

1.43E-05 

2.82E+08 

-7.14E-06 

1.06E+08 

3.57E^06 

1.42E+08 

-3.57E-06 

l,88E+08 

-3.57E-06 

2.82E+08 

1.43E-05 

1.06E+08 

-7.14E-06 

1.42E+08 

-7.14E-06 

1.93E+08 

-7.14E-06 

2.82E+08 

-3.57E-06 

1.08E+08 

5.36E^06 

1.42E+08 

4.43E-05 

1.93E+08 

1.07E-05 

2.82E+08 

O.OOE+OO 

1.08E+08 

7.14E-06 

1.45E+08 

5.36E-06 

1.93E+08 

7.14E-06 

3.04E+08 

3.57E-06 

1.08E+08 

1.07E-05 

1.45E+08 

2.14E-05 

1.93E+08 

-7.14E-06 

3.04E+08 

1.79E-06 

1.08E+08 

1.07E-05 

1.45E+08 

1.43E-05 

1.93E+08 

-3.57E-06 

3.04E+08 

-3.57E-06 

L08E+08 

3.57E-06 

E45E+08 

1.43E-05 

1.93E+08 

3.57E-06 

E18E+08 

-1.07E-05 

1.45E+08 

1.43E-05 

1.93E+08 

-3.57E-06 

1.18E+08 

4.79E-06 

1.53E+08 

-7.14E-06 

2.11E+08 

5.36E-06 
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Hova  SD-2 


Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

2.25E-K)7 

O.OOE+OO 

4.90E+07 

7.14E-06 

6.92E+07 

O.OOE+OO 

9.57E+07 

7.14E-06 

2.25E+07 

O.OOE+00 

4.90E+07 

-7.14E-06 

7.38E+07 

-2.38E-06 

9.57E+07 

7.14E-06 

2.25E+07 

O.OOE+OO 

4.90E+07 

7.14E-06 

7.38E+07 

1.43E-05 

9.57E+07 

1.43E-05 

2.68F,+07 

O.OOE+OO 

4.90E+07 

O.OOE+OO 

7.38E+07 

O.OOE+OO 

1.07E+08 

O.OOE+OO 

2.68E.+07 

O.OOE+OO 

4.90E+07 

O.OOE+OO 

7.38E+07 

O.OOE+OO 

1.07E+08 

O.OOE+OO 

2.68E+07 

7.14E-06 

4.90E+07 

O.OOE+OO 

7.38E+07 

7.14E-06 

1.07E+08 

7.14E-06 

2.68U,+07 

O.OOE+OO 

5.46E+07 

O.OOE+OO 

7.38E+07 

O.OOE+OO 

1.07E+08 

7.14E-06 

2.68E-f07 

7.14E-06 

5.46E+07 

O.OOE+OO 

7.38E+07 

7.14E^06 

1.07E+08 

1.43E-05 

3.17E+07 

3.57E-06 

5.46E+07 

O.OOE+OO 

7.70E+07 

7.14E-06 

1.08E+08 

1.43E^05 

3.17E+07 

7.14E-06 

5.46E+07 

7.14E-06 

7.70E+07 

7.14E-06 

1.08E+()8 

1.43E^05 

3.17E+07 

3.57E-06 

5.46E+07 

O.OOE+OO 

7.70E+07 

O.OOE+OO 

1.08E+08 

1.43E~05 

3.17E+07 

7.14E-06 

5.46E+07 

O.OOE+OO 

8.11E+07 

O.OOE+OO 

1.08E+08 

7.14E-06 

3.17E+07 

3.57E-06 

5.56E+07 

3.57E-06 

8.11E+07 

O.OOE+OO 

1.08E+08 

3.57E-06 

3.17E4-07 

O.OOE+OO 

5.56E+07 

O.OOE+OO 

8.11E+07 

3.57E-06 

1.08E+08 

1.43E-05 

3.70E4()7 

7.14E-06 

5.56E+07 

O.OOE+OO 

8.11E+07 

7.14E-06 

l.lOE+08 

-2.38E  06 

3.70E+()7 

7.14E-06 

5.56E+07 

-7.14E-06 

8.11E+07 

-3.57E-06 

l.lOE+08 

7.14E-06 

3.70E+07 

O.OOE+OO 

S.56E+07 

O.OOE+OO 

8.11E+07 

-7.14E-06 

l.lOE+08 

-7.14E-06 

3.70E+07 

O.OOE+OO 

5.56E+07 

O.OOE+OO 

8.45E+07 

O.OOE+OO 

l.lOE+08 

-3.57E-06 

3.70E+07 

O.OOE+OO 

5.56E+07 

1.43E-05 

8.45E+07 

O.OOE+OO 

l.lOE+08 

3.57E-06 

3.70E-t-07 

O.OOE+OO 

5.75E+07 

O.OOE+OO 

8.45E+07 

O.OOE+OO 

l.lOE+08 

O.OOE+OO 

3.70E4-07 

7.14E-06 

5.75E+07 

O.OOE+OO 

8.45E+07 

1.43E-05 

1.10E+08 

O.OOE+OO 

3.70E+()7 

O.OOE+OO 

5.75E+07 

O.OOE+OO 

8.45E+07 

O.OOE+OO 

1.19E+08 

O.OOE+OO 

3.70E+07 

O.OOE+OO 

6.23E+07 

O.OOE+OO 

8.45E+07  ! 

1.43E-05 

1.19E+08 

O.OOE+OO 

3.87E+()7 

-7.14E-06 

6.23E+07 

3.57E-06 

8.45E+07  ! 

O.OOE+OO 

1.19E+08 

O.OOE+OO 

3.87E+07 

O.OOE+OO 

6.23E+07 

O.OOE+OO 

8.45E+07 

O.OOE+OO 

1.25E+08 

3.57E-06 

3.87E4-07 

O.OOE+OO 

6.23E+07 

O.OOE+OO 

9.18E+07 

7.14E-06 

1.25E+08 

O.OOE+OO 

4.28E+07 

1.19E-06 

6.23E+07 

3.57E-06 

9.18E+07 

1.43E^05 

1.25E+08 

-7.14E-06 

4.28Ef07 

7.14E4)6 

6.38E+07 

O.OOE+OO 

9.18E+07 

1.07E-05 

1.25E+08 

3.57E-06 

4.28E+07 

O.OOE+OO 

6.38E+07 

7.14E-06 

9.18E+07 

3.57E-06 

1.25E+08 

-3.57E-06 

4.28E-f07 

O.OOE+OO 

6.38E+07 

-7.14E-06 

9.18E+07 

1.07E-05 

1.25E+08 

3.57E-06 

4.28E+07 

7.14E  06 

6.38E+07 

O.OOE+OO 

9.47E+07 

3.57E-06 

1.25E+08 

-3.57E-06 

4.28E+07 

O.OOE+OO 

6.38E+07 

O.OOE+OO 

9.47E+07 

3.57E-06 

1.25E+08 

-3.57E-06 

4.28E+07 

O.OOE+OO 

6.38E+07 

7.14E-06 

9.47E+07 

-3.57E-06 

1.27E+08 

1.43E-05 

4.28E+07 

7.14E-06 

6.38E+07 

7.14E-06 

9.47E+07 

-7.14E-06 

1.27E+08 

1.43E-0'5 

4.28E+07 

O.OOE+OO 

6.38E+07 

O.OOE+OO 

9.47E+07 

-3.57E-06 

1.27E+08 

1.07E-05 

4.62E+07 

-7.14E-06 

6.38E+07 

O.OOE+OO 

9.47E+07 

3.57E-06 

1.27E+08 

O.OOE+OO 

4.62E+07 

-7.14E^06 

6.87E+07 

-7.14E-06 

9.47E+07 

7.14E-06 

1.27E+08 

O.OOE+OO 

4.62E+07 

O.OOE+OO 

6.87E+07 

-3.57E-06 

9.47E+07 

-7.14E~06 

1.27E+08 

7.14E^06 

4.62E+07 

-7.14E-06 

6.87E+07 

O.OOE+OO 

9.47E+07 

-3.57E-06 

1.27E+08 

1.43E^05 

4.62E+07 

O.OOE+OO 

6.87E+07 

O.OOE+OO 

9.57E+07 

O.OOE+OO 

1.27E+08 

3.57E-06 

4.90E+07 

7.14E-06 

6.87E+07 

O.OOE+OO 

9.57E+07 

O.OOE+OO 

1.27E+08 

O.OOE+OO 

4.9()E+07 

7.14E-06 

6.92E+07 

3.57E-06 

9.57E+07 

7.14E-06 

1.42E+08 

1.07E-05 

4.90E+07  ! 

O.OOE+OO 

6.92E+07 

7.14E-06 

9.57E+07 

7.14E-06 

1.42E+08 

O.OOE+OO 
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Hoya  SD-2  (continued) 


Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

Dose,  rad 

Strain 

1.42E+08 

O.OOE-fOO 

1.88E+08 

O.OOE+OO 

1.94E+08 

7.14E-06 

2.85E+08 

2.14E-05 

1.42E4-08 

O.OOE+00 

1.88E+08 

O.OOE+OO 

1.94E+08 

7.14E-06 

2.85E+08 

1.07E-05 

1.42E+08 

O.OOE+00 

1.88E+08 

2.14E-05 

1.94E+08 

1.43E-05 

2.85E+08 

1.79E-05 

1.42E+08 

O.OOE+OO 

l,88E+08 

O.OOE+OO 

1.94E+08 

O.OOE+OO 

2.85E+08 

2.50E-05 

1.42E+08 

1.07E-05 

1.88E+08 

O.OOE+OO 

1.94E+08 

O.OOE+OO 

2.85E+08 

1.07E-05 

1.47E+08 

4.76E~06 

l,88E+08 

7.14E-06 

1.94E+08 

1.43E~05 

2.85E+08 

O.OOE+OO 

1.47E+08 

1.43E-05 

1.88E+08 

O.OOE+OO 

2.13E+08 

7.14E-06 

2.85E+08 

3.57E-06 

1.47E+08 

3.57E-06 

1.88E+08 

O.OOE+OO 

2.13E+08 

1.43E-05 

2.85E+08 

7.14E-06 

1.47E+08 

7.14E-06 

1.88E+08 

1.43E-05 

2.13E+08 

1.07E-05 

2.85E+08 

-3.57E-06 

1.47E+08 

1.43E-05 

1.88E+08 

O.OOE+OO 

2.13E+08 

1.43E-05 

2.85E+08 

1.43E-05 

E47E+08 

7.14E-06 

1.88E+08 

7.14E-06 

2.13E+08 

1.43E-05 

2.85E+08 

1.07E-05 

1.47E+08 

7.14E^06 

1.88E+08 

7.14E-06 

2.37E+08 

7.14E-06 

2.85E+08 

7.14E-06 

1.60E+08 

O.OOE+OO 

1.90E+08 

2.14E~05 

2.37E+08 

2.14E-05 

2.88E+08 

1.43E-05 

1.60E+08 

7.14E-06 

1.90E+08 

1.43E-05 

2.37E+08 

1.79E-05 

2.88E+08 

1.43E-05 

E60E+08 

-3.57E-06 

1.90E+08 

7.14E-06 

2.47E+08 

1.43E-05 

2.88E+08 

-3.57E-06 

1.60E+08 

O.OOE+OO 

1.90E+08 

3.57E-06 

2.47E+08 

1.43E-05 

2.88E+08 

3.57E-06 

1.60E+08 

-3.57E-06 

1.90E+08 

7.14E-06 

2.47E+08 

1.55E-05 

2.88E+08 

7.14E-06 

1.67E+08 

-7.14E-06 

1.90E+08 

1.07E-05 

2.60E+08 

1.07E-05 

2.88E+08 

1.43E-05 

1.67E+08 

1.43E-05 

1.90E+08 

2.14E-05 

2.60E+08 

-3.57E-06 

2.88E+08 

7.14E-06 

1.67E+08 

2.62E-05 

1.92E+08 

1.43E-05 

2.60E+08 

O.OOE+OO 

2.88E+08 

1.07E-05 

1.68E+08 

1.43E-05 

1.92E+08 

O.OOE+OO 

2.60E+08 

3.57E-06 

2.88E+08 

2.14E-05 

1.68E+08 

2.14E-05 

1.92E+08 

2.14E-05 

2.60E+08 

1.43E-05 

2.88E+08 

3.57E-06 

1.68E+08 

O.OOE+OO 

1.92E+08 

-7.14E-06 

2.71E+08 

1.07E-05 

2.88E+08 

O.OOE+OO 

1.68E+08 

1.07E-05 

1.92E+08 

2.86E-05 

2.71E+08 

O.OOE+OO 

2.88E+08 

3.57E-06 

1.68E+08 

O.OOE+OO 

1.92E+08 

2.86E-05 

2.71E+08 

-3.57E-06 

2.88E+08 

1.43E  05 

1.68E+08 

1.43E-05 

1.92E+08 

O.OOE+OO 

2.71E+08 

-3.57E-06 

2.88E+08 

1.07E-05 

1.68E+08 

3.57E-06 

1.92E+08 

7.14E-06 

2.71E+08 

3.57E-06 

2.88E+08 

7.14E-06 

1.68E+08 

E43E-05 

1.92E+08 

7.14E-06 

2.71E+08 

7.14E-06 

2.88E+08 

1.43E-05 

1.75E+08 

O.OOE+OO 

1.92E+08 

7.14E-06 

2.71E+08 

1.43E-05 

3.31E+08 

1.43E-05 

1.75E+08 

O.OOE+OO 

1.92E+08 

-7.14E-06 

2.79E+08 

1.43E~05 

3.31  E+08 

2.14E-05 

1.75E+08 

7.14E-06 

1.92E+08 

O.OOE+OO 

2.79E+08 

O.OOE+OO 

3.31E+08 

2.62E-05 

1.75E+08 

7.14E-06 

1.92E+08 

-7.14E-06 

2.79E+08 

1.07E-05 

3.48E+08 

1.43E-05 

1.75E+08 

O.OOE+OO 

1.92E+08 

1.43E-05 

2.79E+08 

-3.57E-06 

3.48E+08 

7.14E-06 

1.77E+08 

3.57E-06 

1.94E+08 

2.14E-05 

2.79E+08 

1.07E-05 

3.48E+08 

3.57E-06 

1.77E+08 

O.OOE+OO 

1.94E+08 

O.OOE+OO 

2.79E+08 

3.57E-06 

3.48E+08 

1.43E-05 

1.77E+08 

1.43E-05 

1.94E+08 

7.14E-06 

2.79E+08 

O.OOE+OO 

3.48E+08 

1.79E-05 

1.82E+08 

1.43E-05 

1.94E+08 

7.14E-06 

2.79E+08 

-3.57E-06 

3.63E+08 

1.43E-05 

1.82E+08 

O.OOE+OO 

l,94E+08 

O.OOE+OO 

2.79E+08 

7.14E-06 

3.63E+08 

7.14E-06 

E82E+08 

O.OOE+OO 

1.94E+08 

2.14E-05 

2.79E+08 

1.43E-05 

3.63E+08 

7.14E-06 

1.82E+08 

O.OOE+OO 

1.94E+08 

1.43E-05 

2.79E+08 

7.14E^06 

3.63E+08 

1.43E-^05 

1.82E+08 

O.OOE+OO 

1.94E+08 

2.14E-05 

2.79E+08 

7.14E-06 

3.63E+08 

7.14E-06 

1.82E+08 

O.OOE+OO 

1.94E+08 

2.86E-05 

2.85E+08 

3.57E-06 

3.63E+08 

1.43E-05 

E82E+08 

2.14E-05 

1.94E+08 

7.14E-06 

2.85E+08 

r  7.14E-06  1  3.63E+08 

2.14E-05 
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Dose,  rad 

Strain 

3.73E+08 

2.14E  05 

3.73E+()8 

2.14E-05 

3.73E+()8 

2.50E-05 

3.74E-i-08 

2.14E-05 

3.74E+08 

1.43E-05 

3.74E+08 

1.07E-05 

3.74E+08 

7.14E-06 

3.74E+08 

7.14E-06 

■ 

3.74E-f08 

2.14Ed)5 

1 

3.74E+08 

2.86E-05 

1 

3.74E+08 

2.14E-05 

1 

3.74E-i'()8 

1.43E-05 

1 

3.81E-H)8 

1.43E-05 

1 

3.81E+08 

2.14E-05 

3.81E+()8 

2.86E4)5 

3.81  E+08 

2.14Ed)5 

3.81E+08 

2.86E-05 

3,81  E+08 

2.86E-05 

3.81  E408 

1.43Edl5 

3.81  E+08 

7.14E-06 

3.81E+08 

1.43E-05 

3.81  E+08 

1.43E  05 

3.81  E+08 

7.14E  06 

3.81E+08 

2.86E-05 

■ 

3.81  E+08 

2.50E-05 

1 

3.81E+0iP 

2.14E-05 

■ 

3.85E+08 

2.14E^05 

3.85E+08 

2.14E-05 

3.85E+08 

1.43E-05 

3.85E+08 

7.14E^06 

3.85E+08 

1.79E-05 

3.85E+08 

2.86E-()5 

3.85E+08 

1.43E-05 

■ 

3.85E+08 

1.43E-05 

I 

3.85E+08 

2.86E-05 

■ 

3.85E+08 

7.14E  06 

3.85E+08 

1.43E-05 

3.85E+08 

1.43E-05 

3.85E+08 

2.14E-05 

3.85E+08 

2.50E-^05 

■ 

3.85E+08 

1.43E-05 

1 

3.85E+08 

2.14E-05 

■ 

APPENDIX  B 

Sink/ float  measurements 


Pyrex 


27  Jun  00 

Sink,  °C 

Ref  sink,  °C 

AT 

a  liquid 

Ad/d 

A1./L 

A 

32.34 

32.29 

0.05 

8.40E-04 

-4.20E-05 

1.40E-05 

B 

32.32 

32.29 

0.03 

8.40E-04 

-2.52E-05 

8.40E-06 

C 

31.72 

.32.26 

-0.54 

8.40E-04 

4.54E-04 

-1.51E-04 

D 

31.13 

32.39 

-1.26 

8.40E-04 

1.06E-03 

-3.53E-04 

E 

28.33 

32.19 

-3.86 

8.40E-04 

3.24E-03 

-1.08E-03 

28  ]un  00 

A 

32.24 

32.19 

0.05 

8.40E-04 

-4.20E-05 

1.40E-05 

B 

32.18 

32.53 

-0.35 

8.40E-04 

2.94E-04 

-9.80E-05 

C 

31.87 

32.04 

-0.17 

8.40E-04 

1.43E-04 

-4.76E-05 

D 

30.84 

32.35 

-1.51 

8.40E-04 

1.27E-03 

-4.23E-04 

E 

28.14 

31.96 

-3.82 

8.40E-04 

3.21  E-03 

-1.07E-03 

29  |un00#l 

A 

32.29 

32.31 

-0.02 

8.40E-04 

1.68E-05 

-5.60E-06 

B 

32.24 

32.24 

0 

8.40E-04 

O.OOE+00 

O.OOE-iOO 

C 

32.02 

32.39 

-0.37 

8.40E-04 

3.11E-04 

-1.04E-04 

D 

30.94 

32.36 

-1.42 

8.40E-04 

1.19E-03 

-3.98E-04 

E 

28.49 

32.35 

-3.86 

8.40E-04 

3.24E-03 

-1.08E-03 

29  )un  00  #2 

A 

32.43 

32.31 

0.12 

8.40E-04 

-l.OlE-04 

3.36E-05 

B 

32.24 

32.36 

-0.12 

8.40E-04 

l.OlE-04 

-3.36E-05 

C 

31.89 

31.98 

-0.09 

8.4aE-04 

7.56E-05 

-2.52E-05 

D 

31 

32.4 

-1.4 

8.40E-04 

1.18E-03 

-3.92E-04 

E 

28.65 

32.12 

-3.47 

8.40E-04 

2.91E-03 

-9.72E-04 

30)un00#l 

A 

32.27 

32.28 

-0.01 

8.40E-04 

8.40E-06 

-2.80E-06 

B 

32.25 

32.02 

0.23 

8.40E-04 

-1.93E-04 

6.44E-05 

C 

31.74 

31.94 

-0.2 

8.40E-04 

1.68E-04 

-5.60E-05 

D 

30.93 

31.89 

-0.96 

8.40E-04 

8.06E-04 

-2.69E-04 

E 

28.26 

32.36 

-4.1 

8.40E-04 

3.44E-03 

-1.15E-03 

30  |un  00  #2 

A 

32.24 

32.12 

0.12 

8.40E-04 

-1.01  E-04 

3.36E-05 

B 

32.3 

32.17 

0.13 

8.40E-04 

-1.09E-04 

3.64E-05 

C 

31.86 

31-95 

-0.09 

8.40E-04 

7.56E-05 

-2.52E-05 

D 

31.03 

31.96 

-0.93 

8.40E-04 

7.81E-04 

-2.60E-04 

E 

28.15 

32.46 

-4.31 

8.40E-04 

3.62E-03 

-1.21E-03 

AVERAGE 

A 

32.30 

32.25 

5.17E-02 

8.40E-04 

-4.34E-05 

1.45E-05 

B 

32.26 

32.27 

-1.33E-02 

8.40E-04 

1.12E-05 

-3.73E-06 

C 

31.85 

32.09 

-2.43E-01 

8,40E-04 

2.04E-04 

-6.81  E-05 

D 

30.98 

32.23 

-1.25E+00 

8.40E-04 

1.05E-03 

-3.49E-04 

E 

28.34 

32.24 

-3.90E+00 

8.40E-04 

3.28E-03 

-1.09E-03 

Standard  Dev 

A 

0.07 

0.08 

6.05E-02 

5.08E-05 

1.69E-05 

B 

0.05 

0.17 

2.03E-01 

1.71  E-04 

5.69E-05 

C 

0.11 

0.19 

1.78E-01 

1.49E-04 

4.98E-05 

D 

0.10 

0.23 

2.47E-01 

2.08E-04 

6.92E-05 

E 

0.20 

0.18 

2.84E-01 

2.38E-04 

7.94E-05 

avg  std  dev 

5.45E-05 

std  dev  of  avg 

2.22E-05 
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Hova  SD-2 


5  )ul00#l 

Sink,  °C 

Ref  sink,  °C 

AT 

a  liquid 

Ad/d 

AI,/E  1 

7\ 

22.87 

22.87 

0 

8.40E-04 

O.OOE+00 

O.OOE-fOO 

B 

22.78 

22.88 

-0.1 

8.40E-04 

8.40E-05 

-2.80E-05 

C 

22.66 

22.76 

-0.1 

8.40E-04 

8.40E-05 

-2.80E-05 

D 

22.78 

22.69 

0.09 

8.40E-04 

-7.56E-05 

2.52E-05 

E 

22.53 

22.79 

-0.26 

8.40E-04 

2.18E-04 

-7.28E-05 

5  |ul  00  #2 

A 

22.77 

22.79 

-0.02 

8.40E-04 

1.68E-05 

-5.60E-06 

B 

22.8 

22.77 

0.03 

8.40E-04 

-2.52E-05 

8.40E-06 

C 

22.8 

22.77 

0.03 

8.40E-04 

-2.52E-05 

8.40E-06 

D 

22.69 

22.74 

-0.05 

8.40E-04 

4.20E-05 

-1.40E-05 

E 

22.52 

22.66 

-0.14 

8.40E-04 

l,18E-04 

-3.92E-05 

5  |ul  00  #3 

A 

22.81 

22.72 

0.09 

8.40E-04 

-7.56E-05 

2.52E-05 

B 

22.91 

22.81 

0.1 

8.40E-04 

-8.40E-05 

2.80E-05 

C 

22.85 

22.88 

-0.03 

8.40E-04 

2.52E-05 

-8.40E-06 

D 

22.8 

22.75 

0.05 

8.40E-04 

-4.20E-05 

1.40E-05 

E 

22.69 

22.81 

-0.12 

8.40E-04 

l.OlE-04 

-3.36E-05 

6  )ul  00  #1 

A 

22.71 

22.74 

-0.03 

8.40E-04 

2.52E-05 

-8.40E-06 

B  1 

22.66 

22.68 

-0.02  n 

8.40E-04 

1.68E-05 

-5.60E-06 

C 

22.71 

22.69 

0.02 

8.40E-04 

-1.68E-05 

5.60E-06 

D 

22.74 

22.69 

0.05 

oo 

o 

-4.20E-05 

1.40E-05 

E 

22.56 

22.66 

-0.1 

8.40E-04 

8.40E-05 

-2.80E-05 

6  lul  00  #2 

A 

22.77 

22.77 

0 

8.40E-04 

O.OOE+00 

O.OOE-iOO 

B 

22.78 

22.73 

0.05 

8.40E-04 

-4.20E-05 

1.40E-05 

C 

22.85 

22.82 

0.03 

8.40E-04 

-2.52E-05 

8.40E-06 

D 

22.74 

22.72 

0.02 

8.40E-04 

-1.68E-05 

5.60E-06 

E 

22.01 

22.76 

-0.15 

8.40E-04 

1.26E-04 

-4.20E-05 

6  |ul  00  #3 

A 

22.71 

22.7 

0.01 

8.40E-04 

-8.40E-06 

2.80E-06 

B 

22.68 

22.6 

0.08 

8.40E-04 

-6.72E-05 

2.24E-05 

C 

22.83 

22.8 

0.03 

8.40E-04 

-2.52E-05 

8.40E-06 

D 

22.68 

22.68 

0 

8.40E-04 

O.OOE-l-OO 

O.OOE+00 

E 

22.61 

22.68 

-0.07 

8.40E-04 

5,88E-05 

-1.96E-05 

AVERAGE 

A 

22.77 

22.77 

8.33E-03 

8.40E-04 

-7.00E-06 

2.33E-06 

B 

22.77 

22.75 

2.33E-02 

8.40E-04 

-1.96E-05 

6.53E-06 

C 

22.78 

22.79 

-3.33E-03 

8.40E-04 

2.80E-06 

-9.33E-07 

D 

22.74 

22.71 

2.67E-02 

8.40E-04 

-2.24E-05 

7.47E-06 

E 

22.59 

22.73 

-1.40E-01 

8.40E-04 

1.18E-04 

-3.92E-05 

Stnncinrcl  Dev 

A 

0.06 

0.06 

4.2613-02 

O.OOE+00 

3.58E-05 

1.19E-05 

B 

0.09 

0.10 

7.34E-02 

O.OOEM-00 

6.17E-05 

2.06E-05 

c: 

0.08 

0.06 

5.28E-02 

O.OOE+00 

4.43E-05 

1.48E-05 

D 

0.05 

0.03 

4.84E-02 

O.OOE-i-OO 

4.07E-05 

1.36E-05 

E 

0.06 

0.07 

6.54E-02 

O.OOE+00 

5.50E-05 

1.83E-05 

std  dev 

1.58E-05 

std  dev  of  avg 

6.46E-06 

100 
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